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Abstract. We present a genealogy for superprocesses with a non-homogeneous quadratic 
branching mechanism, relying on a weighted version of the superprocess and a Girsanov the- 
orem. We then decompose this genealogy with respect to the last individual alive (William's 
decomposition). Letting the extinction time tend to infinity, we get the Q-process by looking at 
the superprocess from the root, and define another process by looking from the top. Examples 
including the multitype Feller diffusion and the superdiffusion are provided. 



1. Introduction 

Even if superprocesses with very general branching mechanisms are known, most of the works 
devoted to the study of their genealogy are concerned with homogeneous branching mechanisms, 
that is, populations with identical individuals. Four distinct approaches have been proposed for 
describing these genealogies. When there is no spatial motion, superprocesses are reduced to 
continuous state branching processes, whose genealogy can be understood by a flow of subordi- 
nators, see Bertoin and Le Gall [5], or by growing discrete trees, see Duquesne and Winkel |13j . 
With a spatial motion, the description of the genealogy can be done using the lookdown process 
of Donnelly and Kurtz [IT] or the snake process of Le Gall [22]. Some works generalize both 
constructions to non- homogeneous branching mechanisms: Kurtz and Rodriguez [20] recently 
extended the lookdown process in this direction whereas Dhersin and Serlet proposed in |10] 
modifications of the snake. 

Using the genealogy, it is natural to consider the corresponding Williams' decomposition, 
which is named after the work of Williams |32] on the Brownian excursion. After Aldous 
recognized in [3] the genealogy of a branching process in this excursion, they also designate de- 
compositions of branching processes with respect to their height, see Serlet [31] for the quadratic 
branching mechanism or Abraham and Delmas [1] for general branching mechanism. Their in- 
terest is twice: they allow to understand the behavior of processes at the top, see Goldschmidt 
and Haas [18] for an application of this approach, and to investigate the process conditioned on 
non extinction, or Q-process, see |31] and Overbeck |25j . 

For Markov processes with absorbing states, the Q-process is defined as the process condi- 
tioned on non absorption in remote time, see Darroch and Seneta [9]. Lamperti and Ney [21] 
found a simple construction in the case of discrete branching processes. Later on, Roelly and 
Rouault [28] provided a superprocess version of this result. Q-processes have intrinsic interest 
as a model of stochastic population, see Chen and Delmas [7]. They also find application in 
the study of the associated martingale, see Lyons, Pemantle and Peres [24| in a discrete setting. 
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Understanding this martingale allows to better understand the original process, see Englander 
and Kyprianou [16] for superprocesses with non homogeneous branching mechanism. 

Our primary interest is to present a genealogy for superprocess with a non-homogeneous 
quadratic branching mechanism, to condition it with respect to its height (this is the William's 
decomposition), and to study the associated Q-process. 

Let X = (Xt,t > 0) be an (£, /3, a) superprocess over a Polish space E. The underlying spatial 
motion Y = ( Yi , t > 0) is a Markov process with infinitesimal generator C started at x under P x . 
The non-homogeneous quadratic branching mechanism is denoted by ip(x, A) = j3{x)\ + a(x)A 2 , 
for suitable functions /3 and a (explicit conditions can be found in Section [2]) . Let be the 
distribution of X started from the finite measure v on E, and N x be the corresponding canonical 
measure of X with initial state x. In particular, the process X under ¥ v is distributed as 
X^gx X % , where Yliel (dX) is a Poisson Point measure with intensity J xeE v(dx)N x (dX). We 
define the extinction time of X: H mSuX = inf{t > 0,Xt = 0}, and assume that X suffers almost 
sure extinction, that is ~N X [H mSuX = oo] = for all x 6 E. Using an /i-transform from Englander 
and Pinsky [17] and a Girsanov transformation from Perkins [22], we provide a genealogical 
structure for the superprocess X, see Proposition 13.12] by transferring the genealogical structure 
of an homogeneous superprocess. 

We define the function Vh(x) = N x [Xh ^ 0] = N x [H max > h] and a family of probability 
measures by setting: 

V o < t < h, dP<x ' Pf = dhVh -^ c -fods a A ^(r.,« h _.(y.)) 
dP x \ Vt d h v h {x) 

where T>t = cr(Y s ,0 < s < t) is the natural filtration of Y, see Lemma 14.101 Using the ge- 
nealogical structure of X, we give a decomposition of the superprocess X with respect to an 
individual chosen at random, also called a Bismut decomposition, in Proposition 14.41 The fol- 
lowing Theorem, see Corollary 14.131 for a precise statement and Corollary 14.141 for a statement 
under F u , gives a Williams' decomposition of X, that is a spine decomposition with respect to 
its extinction time i? max . 

Theorem. (Williams' decomposition under ~N X ) Assume that the (C,f3, a) superdiffusion X 
suffers almost sure extinction and some regularities on a and j3. 

(i) The distribution of H max under N x is characterized by: N x [H max > h] = Vh(x). 

(ii) Conditionally on {H max = ho}, the (C,f3,a) superdiffusion X under N x is distributed 
as X( h °^ constructed as follows. Let x E E and Yiq^^ be distributed according to P x h °\ 
Consider the Poisson point measure M = X^jeJ ^(sj, xj ) on IPj^o) x ^ with intensity: 

21 [o,fc )( s ) ds 1 {H max (x)<h () ~s}a(Y s ) N Ys [dX}. 
The process X^ = (X { t h °\t > 0) is then defined for all t > by: 

x ( t h0) = E x L r 

jeJ, sj<t 

The proof of this Theorem relies on a William's decomposition of the genealogy of X, see 
Theorem 14.121 Notice it also implies the existence of a measurable family (N^ , h > 0) of 
probabilities such that is the distribution of X under N x conditionally on {H max = h}. 

We shall from now on consider the case of Y a diffusion on W K or a pure jump process on a 
finite state space. The generalized eigenvalue Aq of the operator /3 — C is defined in Pinsky [27] 
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for diffusion on K d . For finite state space, it reduces to the Perron Frobenius eigenvalue, see 
Seneta [30]. In both cases, we have: 



A = swp{£ €R,3u£ V(C),u > such that ((3 -C)u = l u}- 

We assume that the space of positive harmonic functions for (/? — Ao) — C is one dimensional, 
generated by a function (pQ. From these assumptions, we have that the space of positive harmonic 
functions of the adjoint of (/3 — Ao) — £ is one dimensional, and we denote by cj>o a generator of 
this space. We also assume that (j>o is bounded from below and above by positive constants and 
that the operator (/5 — Ao) — C is product critical, that is J E dx 4>o(x) 4>o(x) < oo. Thanks to the 
product-critical property, the probability measure P^°, given by: 

v Q dP i°\Vt _ MYt) ft ds (0(Y S )-X O ) 

- ' dP x \ Vt - MYo) 

defines a recurrent Markov process (in the sense given by (|7U|) ). Since 0q is bounded from 
below and from above by two positive constants, the non-negativity of Ao implies the weak 
convergence of the spine that is the weak convergence of towards P^ 00 "* which is given by 
Px° , see Proposition 16.81 An explicit expression is given for P^ and Px° in Lemmas 17.41 and 17.71 
The non- negativity of Ao implies the almost sure extinction of X, see Lemma 16.21 Under very 
general conditions, the weak convergence of the spine implies the convergence of the superprocess 
(Corollary 15. 8|) and its genealogy (Theorem I5.5p . We can easily state them in the particular case 
of an underlying motion being a diffusion or a pure jump process on a finite state space. We 

also see that Ni , defined below, is actually the law of the Q-process, defined as the weak limit 



of the probability measures Ni~^ = N x [ • \H max > h], see also Lemma 15. 11 

Theorem (Q-process under N^). Assume that Ao > 0. Let Y be distributed according to Pt° , 
and, conditionally on Y, let M = Y^jex\sj,xi) be a Poisson point measure with intensity: 

21 R+ (s)ds a(Y s )N Ys [dX]. 

Consider the process 

X {oo) = (X t (oo) ,t > 0), which is defined for all t > by: 



^t — Sj 

Sj<t 



and denote by'N^ its distribution. Then, for allt > 0, the distribution of(X s ,s G [0, t]) under 
or Ni~^ converges weakly to (xi°°\s G [0, t]). 

Notice those results also hold under (see Corollary I5.8P . It is interesting to notice that the 
law of the spine P^ is quite different from that of the backbone given in |17| , see also Remark 



Remark 1.1. As noticed by Li [23] . the multitype Dawson Watanabe superprocess can be under- 
stood as a single non-homogeneous superprocess on an extended space. The above Theorem on 
Q-process provides a construction of the Q-process associated to a multitype Dawson Watanabe 
superprocess considered in Champagnat and Roelly [6], and this construction gives a precise 
meaning to "the interactive immigration" introduced in Remark 2.8 of [Bj. 



Remark 1.2. In |16j . a spinal decomposition of the semi-group of a Doob /i-transform of the 
(£,/3,a) superdiffusion is provided, see Theorem 5 of [16]. The second item of Corollary 15.81 
together with Lemma 15 . 1 1 gives a pathwise decomposition of the genealogy and establishes that 
the process considered actually is the Q-process. 
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We also prove the weak convergence of the probability measures (N x , h > 0) backward from 
the extinction time. Let p( -ft -) denote the push forward probability measure of PW, defined by: 

P { - h \(Y s ,s e [-M]) g .) = P w (ft + ^e [-MD e •). 

The product criticality assumption yields the existence of a probability measure p(~°°) such 
that for all x E E, t > 0, and / bounded measurable: 

s e H, o])] > e(-°°) [/(y s , s e [-t, o])] . 

Once again, the convergence of the spine implies the convergence of the superprocess. The 
following result corresponds to the second item of Theorem 15.91 

Theorem (Asymptotic distribution at the extinction time). Assume that Xq > 0. Then the 
process (Xh+ S ,s E [— 1,0]) under weakly converges towards X_ t ™j, where for s < 0: 

x (-oo )= £ ^ 

jeJ, sj<s 

and conditionally on Y with distribution p( - °°) ; ^ieJ^s-.X?) ^ s a Poisson point measure with 
intensity: 

{ ff ram ( A" ) < — s} 

N Ys [dX]. 

Remark 1.3. Considering a superprocess with homogeneous branching mechanism, the Q-process 
may be easily defined from the well known Q-process for the total mass process (see for instance 
[7] in the case of a general branching mechanism). Thus the recurrence condition imposed on 
the spatial motion is not necessary for Williams decomposition, but it seems more natural in 
order to get the asymptotic distribution at the extinction time. 

Remark 1.4. The genealogy of X defined in Proposition 13. 121 allows us to interpret the following 
probability measure P^ ' as the law of the ancestral lineage of an individual sampled at random 
at height t (see the Bismut decomposition, Proposition I4.4D : 

dP i B \5 t e-Jif*^) 

We prove in Lemma [6.131 that, if cj>o is bounded from below and above by positive constants and 
that the operator (/3 — Ao) — £ is product critical, then the ancestral lineage of an individual 
sampled at random at height t under N x converges as t — > oo to the law of the spine, that is 
P^ 8 '*^ converges weakly to P^ - This Feynman-Kac type penalization result (see Chapter 2 of 
Roynette and Yor [29]) heavily relies on the product criticality assumption, but holds without 
restriction on the sign of Ao • It may be interpreted as an example of the so called globular state 
in random polymers, investigated in Cranston, Koralov and Molchanov [8]. 

Outline. We give some background on superprocesses with a non-homogeneous branching 
mechanism in the Section 2. Section 3 begins with the definition of the /i-transform in the sense 
of Englander and Pinsky, Definition 13.41 goes on with a Girsanov Theorem, Proposition 13.71 
and ends up with the definition of the genealogy, Proposition 13.121 by combining both tools. 
Section 4 is mainly devoted to the proof of the William's decomposition, Theorem 14. 121 By the 
way, we give a decomposition with respect to a randomly chosen individual, also known as a 
Bismut decomposition, in Proposition 14.21 Section 5 gives some applications of the Williams' 
decomposition: We first prove in Lemma 15. II that the limit of the superprocesses conditioned to 
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extinct at a remote time coincide with the Q-process (the superprocess conditioned to extinct 
after a remote time) and actually show in Theorem 15.51 that such a limit exists. We also 
consider in Theorem 15.91 the convergence of the process seen from the top (so, backward from 
the extinction time). All previous results are provided with a set of assumptions. We then give 
in Section [6] sufficient conditions for these assumptions to be valid in term of the generalized 
eigenvector and eigenvalue, then check they hold in Section [7] in two examples: the finite state 
space superprocess (with mass process the multitype Feller diffusion) and the superdiffusion. 

2. Notations and definitions 

This section, based on the lecture notes of Perkins [26], provides us with basic material about 
superprocesses, relying on their characterization via the Log Laplace equation. 
We first introduce some definitions: 

• (E, 5) is a Polish space, B its Borel sigma-field. 

• £ is the set of real valued measurable functions and b£ C £ the subset of bounded 
functions. 

• C(E,M), or simply C, is the set of continuous real valued functions on E, Cb C C the 
subset of continuous bounded functions. 

• D(M + , E), or simply D, is the set of cadlag paths of E equipped with the Skorokhod 
topology, T> is the Borel sigma field on D, and T>t the canonical right continuous filtration 
on D. 

• For each set of functions, the superscript . + will denote the subset of the non-negative 
functions: For instance, b£ + stands for the subset of non negative functions of b£. 

• Aif(E) is the space of finite measures on E. The standard inner product notation will 
be used: for g G £ integrable with respect to M G Aif(E), M{g) = J E M(dx)g(x). 

We can now introduce the two main ingredients which enter in the definition of a superprocess, 
the spatial motion and the branching mechanism: 

• Assume Y = (D, T>, T>t, Yt, P x ) is a Borel strong Markov process. "Borel" means that 
x — > P X (A) is B measurable for all A G B. Let E x denote the expectation operator, and 
(Pt,t > 0) the semi-group defined by: Pt(f)(x) = E x [f(Y t )]. We impose the additional 
assumption that Pt : Cb — > Cb- In particular the process Y has no fixed discontinuities. 
The generator associated to the semi-group will be denoted C Remember / belongs to 
the domain V(j£) of L if / G Cb and for some g G Cb, 

(1) f(Yt) — f( x ) — ds g(Y s ) is a P x martingale for all x in E, 

Jo 

in which case g = C(f). 

• The functions a and j3 being elements of Cb, with a bounded from below by a positive 
constant, the non-homogeneous quadratic branching mechanism ijjP> a is defined by: 

(2) ^' a (x, A) = P(x)\ + a(x)X 2 , 

for all x G E and A G M.. We will just write ip for %jjP' a when there is no possible 
confusion. If a and /3 are constant functions, we will call the branching mechanism (and 
by extension, the corresponding superprocess) homogeneous. 

The mild form of the Log Laplace equation is given by the integral equation, for tf>, f G b£ + , 
t > 0, x G E: 
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(3) u t (x) + E x 



ds ifj(Y s ,u t - s (Y s )) 



o 



f(Yt)+ f ds 
Jo 



Theorem 2.1. ([26\, Theorem II. 5. 11) Let (J),f£ b£ + . There is a unique jointly (in t and x) 
Borel measurable solution u{'^(x) of equation J3J) such that u{'^ is bounded on [0, T] x E for all 
T > 0. Moreover, >0forallt>0. 

We shall write w for u*' when (f> is null. 

We introduce the canonical space of continuous applications from [0,oo) to Aif(E), denoted 
by f2 := C(R + ,Mf(E)), endowed with its Borel sigma field J 7 , and the canonical right continuous 
filtration . Notice that J- = J-^ . 

Theorem 2.2. f|26j. Theorem II. 5. 11) Let u{'^(x) denote the unique jointly Borel measurable 
solution of equation )[3j) such that u{ is bounded on [0,T] x E for all T > 0. There exists a 
unique Markov process X = (Q,,F,Ft,X t , (Fi C,l3 ' a \v G M.f(E))) such that: 

(4) \/<f>, / G b£ + , [e- x *V)-f£*» = e~ v ^ . 
X is called the (£, (3,a)-superprocess. 

We now state the existence theorem of the canonical measures: 

Theorem 2.3. (\2G\ . Theorem II. 7.3) There exists a measurable family of a -finite measures 

(Nx'^ ,a \x G E) on (SI, J 7 ) which satisfies the following properties: If^2j£j0~( x j,xi) * s a Poisson 

point measure on ExQ, with intensity u(dx) ~H x C, ^ ,a \ then ^2j e j X^ is an (£, f3, a)-superprocess 
started at v. 

We will often abuse notation by denoting ¥ u (resp. N x ) instead of p[/''' 3 ' ^ (resp. Ni' C '^' a ' ) ), 
and F x instead of F$ x when starting from 5 X the Dirac mass at point x. 

Let X be a (£., j3, a)-superprocess. The exponential formula for Poisson point measures yields 
the following equality: 

(5) V/ G b£ + , N X0 [1 - ] = - logE^ [e~ x *^ ] = u{(x ), 

where u{ is (uniquely) defined by equation (j4|). 
Denote H max the extinction time of X: 

(6) H max = inf{i > 0; X t = 0}. 

Definition 2.4 (Global extinction). The superprocess X suffers global extinction if ^?v(H max < 
oo) = 1 for all v G Mf(E). 

We will need the the following assumption: 
(HI) The (£, j3, a)-superprocess satisfies the global extinction property. 

We shall be interested in the function 

(7) v t (x)=N x [H mBX >t]. 

We set Voo (x) = lim^oo I vt(x). The global extinction property is easily stated using Voq. 
Lemma 2.5. The global extinction property holds if and only if Voo = 0. 
See also Lemma |4, 91 for other properties of the function v. 
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Proof. The exponential formula for Poisson point measures yields: 

KiHma* <t)= e"^ . 
To conclude, let t goes to infinity in the previous equality to get: 

P„(# max < oo) = e-^~) . 

□ 

For homogeneous superprocesses (a and (3 constant), the function v is easy to compute and 
the global extinction holds if and only /3 is non-negative. Then, using stochastic domination 
argument, one get that a (£, /?, a)-superprocess, with f3 non-negative, exhibits global extinction 
(see [16] p. 80 for details). 



3. A GENEALOGY FOR THE SPATIALLY DEPENDENT SUPERPROCESS 

We first recall (Section 13. lj) the /i-transform for superprocess introduced in |17| and then 
(Section 13. 2p a Girsanov theorem previously introduced in |26j for interactive superprocesses. 
Those two transformations allow us to give a Radon-Nikodym derivative of the distribution of 
a superprocess with non-homogeneous branching mechanism with respect to the distribution of 
a superprocess with an homogeneous branching mechanism. The genealogy of the superprocess 
with an homogeneous branching mechanism can be described using a Brownian snake, see |12j . 
Then, in Section 13.31 we use the Radon-Nikodym derivative to transport this genealogy and get 
a genealogy for the superprocess with non-homogeneous branching mechanism. 

3.1. /i-transform for superprocesses. We first introduce a new probability measure on (D, T>) 
using the next Lemma. 

Lemma 3.1. Let g be a positive function ofD(C) such that g is bounded from below by a positive 
constant. Then, the process ( 9 g^ e ~ ^° ds ^ 9 ^ 9 ^ Ys \t > 0) is a positive martingale under P x . 

We set V g (C) = {v e C b , gv G £>(£)}. 

Proof. Let g be as in Lemma [37TI and / 6 V g (C). The process: 

(f9)(Y t ) - (fg)(x) - J ds C(fg)(Y s ), t > o\ 
is a P x martingale by definition of the generator C Thus, the process: 



V g( x ) Jo g(% 

is a P x martingale. We set: 
(8) M{>° = e" fo * <£a/9){Y.) _ f[x) 

1 g{x) 



t ds ^m,t>o 



1 ds e -J°dr (Cg/ 9 )(Y r ) \£(f9)(Y.) C{g){Y s ) (fg)(Y s 



g{x) g(Y s ) g(x) 

Ito's lemma then yields that the process (Af/' 9 , t > 0) is another V x martingale. Take / constant 
equal to 1 to get the result. □ 
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Let P| denote the probability measure on (D,T>) defined by: 

(9) V t > 0, |Pt = iQQ e" /o rfs . 

dPx \ Vt g{x) 

Note that in the case where g is harmonic for the linear operator C (that is Cg = 0), the 
probability distribution P 9 is the usual Doob /i-transform of P for h = g. 

We also introduce the generator C 9 of the canonical process Y under P 9 and the expectation 
operator E 9 associated to P 9 . 

Lemma 3.2. Let g be a positive function ofD(£) such that g is bounded from below by a positive 
constant. Then, we have T> g (C) C V(C 9 ) and 

VueV 9 (C), C 9 (u) = C{9u) - C{9)u . 
Proof. As, for / G T> g (C), the process 

(M/ ,9 ,i > 0) defined by © is a martingale under P x , we 

get that the process: 



9(5-.: 

is a P| martingale. This gives the result. □ 

Remark 3.3. Let ((t,x) — > g(t,x)) be a function bounded from below by a positive constant, 
differentiable in t, such that .) G f (>C) for each t and ((t, x) — >■ dtg(t,x)) is bounded from 
above. By considering the process (t,Yt) instead of Yt, we have the immediate counterpart 
of Lemma 13.11 for time dependent function g(t,.). In particular, we may define the following 
probability measure on (D,T>) (still denoted V 9 X by a small abuse of notations): 

(10) vt>o, d IhRi = ^ e -f^s^ { s, Ys) 
1 ; dP x[Di g(0,x) 

where C acts on g as a function of x. 

We now define the /i-transform for superprocesses, as introduced in [T7] (notice this does not 
correspond to the Doob /i-transform for superprocesses). 

Definition 3.4. Let X = (Xt,t > 0) be an (C,/3,a) superprocess. For g G b£ + , we define the 
h-transform of X (with h = g) as X 9 = (Xf,t > 0) the measure valued process given for all 
t > by: 

(11) Xf(dx) = g{x)X t {dx). 

Note that (jlip holds point-wise, and that the law of the /i-transform of a superprocess may 
be singular with respect to the law of the initial superprocess. 

We first give an easy generalization of a result in section 2 of j!7j for a general spatial motion. 



Proposition 3.5. Let g be a positive function ofT>(C) such that g is bounded from below by a 
positive constant. Then the process X 9 is a (^C 9 , (~^+^) g ; ag^j -superprocess. 

Proof. The Markov property of X 9 is clear. We compute, for / G b£ + : 
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where, by Theorem 12. 2| u satisfies: 
(12) u t (x)+E x 
which can also be written: 



dr ifj(Y r ,u t -r(Y r )) 



ut{x) + E a 



dr ip(Y r ,u t - r (Xr)) 



dr i/>(Y r ,ut-r(Y r )) 



E x [(fg){Y t )]. 



But (I12p written at time t — s gives: 



u t - s (x) + E X 



t-s 



dr ip(Y r ,u t - s -r(Y r )) 



By comparing the two previous equations, we get: 



u t (x) + E a 



dr i/j(Y r ,ut- r (Y r )) 



E x [(fg)(Y t . s )]. 



Ex [u t - s (Y s [ 



and the Markov property now implies that the process: 

ut- a (Y s )- dr ip(Y r ,Ut- r (Y r )) 
Jo 

with s £ [0, t] is a P x martingale. Ito's lemma now yields that the process: 

U t - s (Y s )e-foM£9/ 9 )(Yr) _ f S dr e -f r du (Cg/ 9 )(Y U ) ^ry riUt _ r{ Y r)) _ (jCg/g)(Y r ) U t . r {Y r )) 



with s £ [0, t] is another P x martingale (the integrability comes from the assumption Cg £ Cb 
and 1/g £ Cb)- Taking expectations at time s = and at time s = t, we have: 



u t (x) + E X 



ds e -Iodr(C g /g)(Y r ) ^ {Ys , Ut _ s (Y s )) - {Cg / g){Y s )u t . s {Y s )) 



E, 



■f dr(C 9/g ){Y r ) Ug){ y t) 



We divide both sides by g{x) and expand tjj according to its definition: 



(-)(*) +E a 
9 



, g{Y s ) 

ds — t^t e 



L jo 



. KMCg/gWr) ( { ag){Y s ){^)\Y s ) + {P--){Y s ){^)(Y s ) 
g{x) V 9 9 9 



E, 



g( y *) f*dr(Cg/g)(Y r ) 



9{x) 



By definition of from @, we get that: 

r-t 



ds ((ag)(Y s )(^) 2 (Y s ) + 09 - ^)(r s )(^)(F s ) 
o V S 5 . 



Eg [/(**)]. 



We conclude from Theorem 12.21 that X 9 is a - c ~^^ 9 } ag)-superprocess. □ 

In order to perform the /i-transform of interest, we shall consider the following assumption. 

(H2) 1/a belongs to £>(£). 

Notice that (H2) implies that aC{l/a) £ Cb- Proposition 13.51 and Lemma [3Tl then yield the 
following Corollary. 
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Corollary 3.6. Let X be an (£, f3,a)-superprocess. Assume (H2). The process X l / a is an 
(£., /3,l)-superprocess with: 

(13) C = C 1/a and p = p - a£(l/a). 

Moreover, for all t > 0, the law P x of the process Y with generator C is absolutely continuous 
on T>t with respect to P x and its Radon-Nikodym derivative is given by: 

K ' dP x]Vt a(Y t ) 

We will note P for the law of X l l a on the canonical space (that is P = P^'^' 1 )) and N for its 
canonical measure. Observe that the branching mechanism of X under P, which we shall write 
ip, is given by: 

(15) tp(x, A) = f3(x) A + A 2 , 

and the quadratic coefficient is no more dependent on x. Notice that ¥ au (X G ■) = ¥' u (aX £ •). 
This implies the following relationship on the canonical measures (use Theorem 12.31 to check it): 

(16) a(x)N x [X G •] =N x [aX G •]. 

Recall that v t (x) = N x [H max > t] = N x [X t / 0]. We set v t (x) = N x [X t / 0]. As a is positive, 
equality (|16p implies in particular that, for all £ > and x £ E: 



(17) a(x)v t (x) = v t (x). 

3.2. A Girsanov type theorem. The following assumption will be used to perform the Gir- 
sanov change of measure. 

(H3) Assume (H2) holds. The function j3 defined in (I13|) is in V(C), with C defined 
in (EES). 



For z S R, we set z + = max(z,0). Under (H2) and (H3), we define: 

(18) & = sup max (p(x),J(0*(x) - 2C0)(x))+) and q(x) = & • 
Notice that q > 0. 

We shall consider the distribution of the homogeneous (C,/3q, l)-superprocess, which we will 
denote by P° (P° = P^-*' 1 )) an d its canonical measure N°. Note that the branching mechanism 
of X under P° is homogeneous (the branching mechanism does not depend on x). We set tp° for 
Since ip° does not depend anymore on x we shall also write ijj W f° r V ;0 ( x ; -^) : 

(19) V°(A) = /3 A + A 2 . 

Proposition 13.71 below is a Girsanov's type theorem which allows us to finally reduce the 
distribution P to the homogeneous distribution P°. We introduce the process M = (M t ,t > 0) 
defined by: 

(20) M t = exp (X (q) - X t (q) - ds X s (if) 
where the function cp is defined by: 

(21) <p(x) = $(x, q{x)) - C(q)(x), x G E. 
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Proposition 3.7. A Girsanov's type theorem. Assume (H2) and (H3) hold. Let X be a 
(£, (3, l)-superprocess. 

(i) The process M is a bounded J- '-martingale under ¥ u which converges a.s. to 

(ii) We have: 

dF° 

—~- = Moo. 

cffiV 

(iii) If moreover (HI) holds, then P®-a.s. we have > 0, the probability measure ¥ u is 
absolutely continuous with respect to P° on T: 

dP?~M> ana dm~ 



We also have: 



(22) q{x) = N° x 



1 



The two first points are a particular case of Theorem IV. 1.6 p. 252 in [26] on interactive drift. 
For the sake of completeness, we give a proof based on the mild form of the Log Laplace equation 
([3]) introduced in Section [H Notice that: 

(23) V°(A) = 4>(x, A + q(x)) - $(x, q(x)). 

Thus, Proposition 13.71 appears as a non-homogeneous generalization of Corollary 4.4 in [2]. We 
first give an elementary Lemma. 

Lemma 3.8. Assume (H2) and (H3) hold. The function ip defined by \21\) is non-negative. 
Proof. The following computation: 

cp(x) = tp(x, q(x)) - C(q){x) = q(x) 2 + j3q(x) - C{q){x) 

^-^y + ~ m ^-f^ -C(g)(x) 



2 



and the definition (|18j) of /3q ensure that the function <p is non-negative. 



□ 



Proof of Proposition \3. 71 First observe that M is J-"-adapted. As the function q also is non- 
negative, we deduce from Lemma 13.81 that the process M is bounded by e x °^ q \ 
Let / G b£ + . On the one hand, we have: 

E x [M t e- Xt{f) ] = t x [e q{x) - Xt{q+f) -ti ds X °M] = ^ 

where, according to Theorem 12.21 rt(x) is bounded on [0,T] x E for all T > and satisfies: 



(24) rt(x) + E, 



ds i;(Y t _ s ,r s (Y t _ s )) 
= E 



ds $(Y t _ s ,q(Y t _ s )) - C(q)(Y t . a )) + (q + f)(Y t ) 
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On the other hand, we have 

E 0[ e -X t (/)j = e -w t {x)^ 

where wt{x) is bounded on [0, T] x E for all T > and satisfies: 

r-t 



w t (x) + E a 



ds ^(Y t _ s ,w s (Y t . s )) 



Using (|23p . rewrite the previous equation under the form: 



ds i>(Y t - a ,(w a + q)(Y t - a )) 



(25) w t (x)+B a 



We now make use of the Dynkin's formula with (H3): 



— E x 


[L 







V x [f(Yt)}. 



ds ?jj(Y t _ s ,q(Y t _ s ))+f(Y t ) 



(26) 



C{q){Y s 



+ E x [q(Yt)}, 



and sum the equations ([25]) and ([26]) term by term to get: 



(27) (w t + q)(x)+E x 



ds i>(Y t - s ,(w s + q)(Y t - s )) 



E, 



ds ^(Y t . s ,q(Y t . s )) - C(q)(Y t . s )) + (q + f)(Y t ) 



The functions r^(x) and Wt(x) + g(x) are bounded on [0, T] x E for all T > and satisfy the 
same equation, see equations (|24]) and (|27|l . By uniqueness, see Theorem l2,14 we finally get that 
Wt + q = r t . This gives: 

(28) E x [M t e- Xt{f) ] = E° x [e- Xt(f) }. 

The Poissonian decomposition of the superprocesses, see Theorem 12.3] and the exponential 
formula enable us to extend this relation to arbitrary initial measures v: 

(29) E u [M t e~ Xt ^] = E°[e- Xt ^]. 

This equality with / = 1 and the Markov property of X proves the first part of item (i). 

Now, a direct induction based on the Markov property yields that, for all positive integer n, 
and fi, . . . , f n e b£ + , < s x < . . . < s n < t: 

(30) Ej,[M t e" El <»<™ Xs * (/l) ] = E°[e _El <'<™ Xs ' (/l) ]. 

And we conclude with an application of the monotone class theorem that, for all non- negative 
J-"t-measurable random variable Z: 

E v [M t Z]=El[Z}. 

The martingale M is bounded and thus converges a.s. to a limit M^. We deduce that for all 
non-negative J-"t-measurable random variable Z: 

(31) E U [M 00 Z]=E° U [Z}. 

This also holds for any non-negative J 7 ^ -measurable random variable Z. This gives the second 
item (ii). 

On {H max < +oo}, then clearly M t converges to e x °W-fo°° ds Xa ^ . Notice that F°(i? ma x = 
+oo) = 0. We deduce from (l3T|) with Z = l{H max =+ao} that F u -a.s. on {H max = +00} , M m = 0. 
This gives the last part of item (i). 
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Now, we prove the third item (iii). Notice that (|3ip implies that P^-a.s. M m > 0. Thanks 
to {HI), we also have that P^-a.s. > 0. Let Z be a non-negative J-^ -measurable random 
variable. Applying ([31]) with Z replaced by 1{jw oo>0 }^/M oo , we get: 



E„[Z] = E v 



Moo 1{M 00 >0}Jj~ 



E, 



L {M oo >0} 



= E° 









This gives the first part of item (iii) . 

Notice that for all positive integer n, and f±, . . . , f n G < s\ < 



< s n , we have 



N 



log (E, 
log (E° 
1 



g5Zl<i<n (/«) 

eEl<i<n-MA)+Jo"~ 



+ 9 



Taking fi = for all i gives (|22|) . This implies: 



]_ _ g5Zl<i< n (fi) 



N 



j +o ° A _ e Ei<i<„^(/0 



The monotone class theorem gives then the last part of item (iii). 



□ 

given by 



3.3. Genealogy for superprocesses. We now recall the genealogy of X under 
the Brownian snake from |12J . We assume (H2) and (H3) hold. 

Let W denote the set of all cadlag killed paths in E. An element w E W is a cadlag path: 
w : [0, T](w)) i?, with 77(10) the lifetime of the path w. By convention the trivial path {x}, 
with x £ E, is a, killed path with lifetime and it belongs to W. The space W is Polish for the 
distance: 

rri(w)Ari(w') 



d( w , w') = 5(w(0),w(0)') + \r)(w) - r](w')\ + 



ds d s {w [0tS] ,w' [QtS] ), 



where d s refers to the Skorokhod metric on the space D([0, s], E), and wj is the restriction of 
w on the interval /. Denote W x the set of stopped paths w such that w(0) = x. We work on 
the canonical space of continuous applications from [0,oo) to W, denoted by f2 := C(M + ,W), 
endowed with the Borel sigma field Q for the distance d, and the canonical right continuous 
filtration Q t = a{W s ,s < t}, where (W s ,s € M + ) is the canonical coordinate process. Notice 
Q = Goo by construction. We set H s = rj(W s ) the lifetime of W s . 



Definition 3.9 (Proposition 4.1.1 and Theorem 4.1.2 of [12]). Fix Wo £ W x - There exists a 
unique W x -valued Markov process W = (CI, Q, Qt, Wt, Pv^ ); called the Brownian snake, starting 
at Wo and satisfying the two properties: 

(i) The lifetime process H = (H s ,s > 0) is a reflecting Brownian motion with non-positive 
drift —Po> starting from Hq = 77(Wo). 

(ii) Conditionally given the lifetime process H, the process (W s ,s > 0) is distributed as 
an inhomogeneous Markov process, with transition kernel specified by the two following 
prescriptions, for < s < s' : 

- W s /(t) = W s (t) for all t < H[ atS t], with F[ s ,y] = inf s < r < s / H r . 

— Conditionally on W a (.ffr a)S /]— ), the path (W s '(Ht s s n + t),0 < t < H s i — Hr s s n) is 
independent of W s and is distributed as Y[o,h s ,-h, s s ,-,) under Pw s (H [s „/]-)■ 

This process will be called the Po-snake started at Wo, and its law denoted by P^ - 



Let ^j^jbixi wi) be a Poisson point measure on E x Q with intensity v(dx) N°[c£W]. 
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We will just write for the law of the snake started at the trivial path {x}. The corre- 
sponding excursion measure of W is given as follows: the lifetime process H is distributed 
according to the Ito measure of the positive excursion of a reflecting Brownian motion with 
non-positive drift — /3o, and conditionally given the lifetime process H, the process (W s ,s > 0) 
is distributed according to (ii) of Definition 13.91 Let 

a = inf{s > 0; H s = 0} 

denote the length of the excursion under N x . 

Let (lg,r > 0, s > 0) be the bicontinuous version of the local time process of H; where l r s 
refers to the local time at level r at time s. We also set w = w(r](w)—) for the left end position 
of the path w. We consider the measure valued process X(W) = (Xt(W),t > 0) defined under 
N£by: 

(32) X t (W)(dx)= r dJt5 Ws (dx). 

Jo 

The /3o-snake gives the genealogy of the (C,/3o, 1) superprocess in the following sense. 

Proposition 3.10 ([12], Theorem 4.2.1). We have: 

The process X(W) is under N° distributed as X under N°. 

Then ^2j e j X(W 3 ) is an (£, /3q, I) -superprocess started at v. 

Notice that, under N^, the extinction time of X(W) is defined by 

mi{t;X t (W) = 0} = sup H s , 

and we shall write this quantity H m . AX or H max (W) if we need to stress the dependence in W. 
This notation is coherent with (J6j). 

We now transport the genealogy of X under N° to a genealogy of X under N. In order to 
simplify notations, we shall write X for X(W) when there is no confusion. 

Definition 3.11. Under (H1)-(H3), we define a measure N x on (£i,G) by: 

Notice the second equality in the previous definition is the third item of Proposition 13.71 
At this point, the genealogy defined for X under will give the genealogy of X under N up 
to a weight. We set 

(33) N x = -L N x . 

a{x) 

Proposition 3.12. We have: 

(i) X under N x is distributed as X under N x . 

(ii) The weighted process X we ^ ht = (x™ eigh \ t > 0) with 

(34) Xr 9h \dx) = [° dJt a(W s )5w (dx), t > 0, 



is under N x distributed as X under N^. 
We may write X wei&ht (W) for x wcl&ht to emphasize the dependence in the snake W. 
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Proof. This is a direct consequence of Definition 13.111 and (|16p . 

We shall say that W under N x provides through (|34p a genealogy for X under 



□ 



4. A Williams' decomposition 

In Section 14,11 we give a decomposition of the genealogy of the superprocesses (£, /3, a) and 
(£, f3, 1) with respect to a randomly chosen individual. In Section 14.21 we & ve a Williams' 
decomposition of the genealogy of the superprocesses (C, (3, a) and (£, /3, 1) with respect to the 
last individual alive. 

4.1. Bismut's decomposition. A decomposition of the genealogy of the homogeneous super- 
process with respect to a randomly chosen individual is well known in the homogeneous case, 
even for a general branching mechanism (see lemmas 4.2.5 and 4.6.1 in |12j). 

We now explain how to decompose the snake process under the excursion measure (N x or N°) 
with respect to its value at a given time. Recall a = inf {s > 0, H s = 0} denote the length of the 
excursion. Fix a real number t S [0, a]. We consider the process Hk» (on the left of t) defined 



on [0,t] by H 
r(s) 



(.<;) 



- -fft for all s G [0, t]. The excursion intervals above of the process 
(.ffi 9 ^ — info< s '< s iJ^f \ < s < t) are denoted Uj g j(9)( c j>4/)- We also consider the process 



Hp 
r {9) 



(on the right of t) defined on [0, a — t) by H, 



id) 



H 



t+s 



Hi. The excursion intervals above 



of the process (Hs — inf < s '< s H^ 1 ' , < s < a — t) are denoted \Jj & j(d)(cj,dj). We define the 
level of the excursion j as sj = Hf- C . if j € jw and Sj = Ht +Cj if j E j( d \ We also define for 
the excursion j the corresponding excursion of the snake: W 3 = (Wi , s > 0) as 

W7'(.) = Wi_ (<y+s)Ad .(. + Si ) ifjGJ^,and Wi(.)=W t+(c . +s)Ad .(. + Si ) ifjGjW. 

We consider the following two point measures on M + x Cl: for e £ {5, d}, 



r(d) 



(35) 



Notice that under (and under N x if (HI) holds), the process W can be reconstructed from 
the triplet (Wt, Rf, Rf). We are interested in the probabilistic structure of this triplet, when t 
is chosen according to the Lebesgue measure on the excursion time interval of the snake. Under 
N°, this result is as a consequence of Lemmas 4.2.4 and 4.2.5 from [12]. We recall this result in 
the next Proposition. 

For a point measure R = J2jej ^(sj,xj) 011 a space MxA? and A C M, we shall consider the 
restriction of R to A x X given by Ra = ^-A(sj)5^ s .^ Xj y 

Proposition 4.1 (|12j. Lemmas 4.2.4 and 4.2.5). For every measurable non-negative function 
F, the following formulas hold: 



(36) 
(37) 



ds F(W s ,R 9 s ,R a s 



e -Ar dr g 



F(Y { 



[0,r)' 



N 



d s l\ F(W S ,R 9 S ,R 



e~ Pot E x 



0,r-)' rL [0,r) 
B,d 



F ( Y [0,t),R[ %i R [O t 



t > 0, 



where under E x and conditionally on Y, R B,g and R B,d are two independent Poisson point 
measures with intensity v B (ds,dW) = ds Ny s [dl4 7 ]. 



The next Proposition gives a similar result in the non-homogeneous case. 
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Proposition 4.2. Under (H1)-(H3), for every measurable non-negative function F , the two 
formulas hold: 



(38) N x 



ds F(W s> R 9 ,R d ) 



dr E* 



,B,d 
'[0.0 



[0,r) 



where under E x and conditionally on Y , Rm^.\ an d ^[oV) are independent Poisson point 
measures with intensity 

(39) z; B (cis, dWQ = (is Ny s [dW] = ds a(Y s )N Ya [dW]; 

and 



(40) N a 



ds a(W s )F(W s ,R 9 s ,R1 



dr Fj x 



[0,r)' -"-[0,r)^ 



where under ~E X and conditionally on Y , R^'f\ an d R^'^ are ^ wo independent Poisson point 



measures with intensity v B . 



[0,0 



Observe there is a weight a(W s ) in (|40p (see also (|34p where this weight appears) which 
modifies the law of the individual picked at random, changing the modified diffusion P x in (|38p 
into the original one P^. 

We shall use the following elementary Lemma on Poisson point measure. 

Lemma 4.3. Let R be a Poisson point measure on a Polish space with intensity v. Let f be 
a non-negative measurable function f such that v[e^ —1) < +oo. Then for any non-negative 
measurable function F , we have: 



(41) 



E 



F(R)e 



R(f) 



E 



F(R) 



'(ef -1) 



where R is a Poisson point measure with intensity u(dx) = v(dx). 

Proof of Proposition \4-%\ We keep notations introduced in Propositions 14,11 and 14.21 We have: 



ds F{W S ,R 9 ,R° S 



N 



e / +o ° ds x s (?) r ds F ( Ws ]R g jR c 

Jo 

ds F(W s ,R 9 ,R d )e^+ Rd ^ 



~P° r dr E 7 



dr E, 



-Por 



dr E 7 



r \ x [Q,r)i -"-[0,0' [0,0'' 

F(Y„ , R B > 9 R B ' d ) e 2 & ds N ^ [ e/o+ °° Xr(W)M - 1 ] 
F(V, ^ K B > 9 U B ' d \ P . 2 fo ds ^) 



dr E 7 







L [0,r)' [0,r) 

where the first equality comes from (HI) and item (iii) of Proposition 13.71 we set f(s,W) 
So 



Jq +0 ° X r (W)(ip) for the second equality, we use Proposition 14.11 for the third equality, we use 
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Lemma [4,31 for the fourth, we use (|22[) for the fifth, and the definition ()18|) of q in the last. This 
proves ([5H]). 

Then replace F{W s ,R 9 s ,R d s ) by a(W s )F(W s ,R 9 s ,R d ) in flM} and use ([11]) as well as d33j) to 
get (SD1). □ 



The proof of the following Proposition is similar to the proof of Proposition 14.21 and is not 
reproduced here. 

Proposition 4.4. Under (H1)-(H3), for every measurable non-negative function F, the two 
formulas hold: for fixed t > 0, 



(42) 



N, 



dJl F(W s ,R 9 ,R d ) 



B,d > 
[0,t)' Jl [0,t)^ 



where under ~E X and conditionally on Y, R B,g and R B,d are two independent Poisson point 
measures with intensity v B defined in 113 9\) . and 



(43) 



N, 



dJt a(W s )F(W s ,R 9 ,R d ) 



E, 



o~fo ds P(Xs) J7(V r , R B >9 r>B,d\ 



where under E x and conditionally on Y, R B ' 9 and R Bid are two independent Poisson point 
measures with intensity v B . 

As an example of application of this Proposition, we can recover easily the following well 
known result. 

Corollary 4.5. Under (H1)-(H3), for every measurable non-negative functions f and g on E, 
we have: 



N, 



X t (f)e 



-x t (g) 



E, 



■Jl da 8^{Y a , Ny,[l-e**-W]) 



In particular, we recover the so-called "many-to-one" formula (with g = in Corollary I4.5P : 



(44) 



N x [X t (f)] = E a 



■fid. P(Y S ) f{Yt) 



Remark 4.6. Equation (j44H justifies the introduction of the following family of probability mea- 
sures indexed by t > 0: 



(45) 



dP 



(B,t) 
x \V t 



■ ;*«fe/9(Y;) 



dP 



x \v t E x 



e -f*ds y9(y.) 



which can be understood as the law of the ancestral lineage of an individual sampled at random 
at height t under the excursion measure N x , and also correspond to Feynman Kac penalization of 
the original spatial motion P x (see [29] ) . Notice that this law does not depend on the parameter 
a. These probability measures are not compatible as t varies but will be shown in Lemma 16.131 
to converge as t — > oo in restriction to P s , s fixed, s < t, under some ergodic assumption (see 
(H9) in Section E]). 



Proof. We set for w £ W with r](w) = t and r±, r2 two point measures on IR + x Q 



F(w, ri ,r 2 ) = f(w)e h ^ +h ^\ 
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where M£ i6 j«W*)) = E Si <t X^\W%. Si (g). We have: 



N, 



X t (f)e 



TVT 

= 


.JO 






= E x 






= Ea; 


e -;*cfa /8(y.) /(y t ) e -/o 2Q ( y J N nN 


x wcight (s) 


= E x 







where we used item (ii) of Proposition 13. 121 for the first and last equality, (|43p with i 7 previously 
defined for the second, formula for exponentials of Poisson point measure and ()33|) for the 
third. □ 



4.2. Williams' decomposition. We first recall the Williams' decomposition for the Brownian 
snake (see [32] for Brownian excursions, |31] for Brownian snake or pQ for general homogeneous 
branching mechanism without spatial motion). 

Under the excursion measures N°, N x and N x , recall that H m£kX = supr 0)CT .i H s . Because of 
the continuity of H, we can define T m ax 

= inf{s > 0, H s = i^max}- Notice the properties of the 
Brownian excursions implies that a.e. H s = -ff max only if s = T max . We set v®(x) = N°[iJ max > t] 
and recall this function does not depend on x. Thus, we shall write v® for v®(x). Standard 
computations give: 

o = A) 

* eA* -1 

The next result is a straightforward adaptation from Theorem 3.3 of pQ and gives the distribution 
of (H max ,W T ^,R°R d T ) under N°. 



Proposition 4.7 (Williams' decomposition under N x ). VFe /iaue: 

(i) T/ie distribution of H max under is characterized by: N x [H max > /i] = v®. 

(ii) Conditionally on {H rnax = ho}, the law ofWT max under N°_ is distributed as V[o,/i ) under 

P*. 

(iii) Conditionally on {H max = ho} and WT max , R? and R T max are under ^x independent 
Poisson point measures on M + x f2 loii/i intensity: 

l [0M) ( S )ds l{ Hmax (W)<h -s} ™W Tma Js) [ dW l 

In other words, for any non-negative measurable function F, we have 



F(H max , WT max , R 9 T „ in ^ R T„, 



d h v° h dh E x 



F(h,Y l0 , h) ,R w ^,R w ^ d ) 



where under E x and conditionally on Y\o,h)i R W '^' 9 and R w >W' d are two independent Poisson 
point measures with intensity i> w ^ h \ds,dW) = lr 0) / l )(s)ds l{H max (W)<h-s} Ny [dW]. 

Notice that items (ii) and (iii) in the previous Proposition implies the existence of a measurable 



family (N I M ft '\ h > 0) of probabilities on (fi, Q) such that N^' w is the distribution of W (more 



o,(fc) 



precisely of (lFr mK , R 9 T ,Rj,)) under N , conditionally on {H n 



h}. 
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Remark 4.8. In Klebaner & al [19], the Esty time reversal "is obtained by conditioning a [dis- 
crete time] Galton Watson process in negative time upon entering state (extinction) at time 
when starting at state 1 at time — n and letting n tend to infinity". The authors then observe 
that in the linear fractional case (modified geometric offspring distribution) the Esty time re- 
versal has the law of the same Galton Watson process conditioned on non extinction. Notice 
that in our continuous setting, the process (H s ,0 < s < T max ) is under a Bessel process 

up to its first hitting time of h, and thus is reversible: (H s ,0 < s < T max ) under N° is dis- 
tributed as (h — flr max -s,0 < s < T max ) under N° . It is also well known (see Corollary 3.1.6 
of [12]) that (H a _ s ,0 < s < a — T max ) under N^'^ is distributed as (H s ,0 < s < T max ) under 
Nx . We deduce from these two points that (X S (1),0 < s < h) under N°' is distributed 
as (Xh_ s (l),0 < s < h) under . This result, which holds at fixed h, gives a pre- limiting 

version of the Esty time reversal in continuous time. Passing to the limit as h — > oo, see Section 
15. 2\ we get the equivalent of the Esty time reversal in a continuous setting. 

Before stating the Williams' decomposition, Theorem I4.12( let us prove some properties for 
the functions vt(x) = N x [i/ max > t] = N x [Xt ^ 0] and vt(x) = Na;[.£f max > t] which will play a 
significant role in the next Section. Recall ()17p states that 



av t = v t . 



Notice also that f)18|) implies that q is bounded from above 
Lemma 4.9. Assume (H1)-(H3). We have: 



by (A) + ID/2. 



(46) q(x) +v\ ?> v t (x) > v^. 

Furthermore for fixed x G E, Vt(x) is of class C 1 in t and we have: 




where the function £ defined by: 



(48) E t (x) = 2{v° t + q(x) - v t (x)) = d x ^{v° t ) - d x i>(x,v t (x)) 



satisfies: 




v t (x) = ® x [X t ^0]=N° X \l {Xtm e/o + °° * > n° x [X t t 0] = vl 



We also have 




{x t =o} 



= q(x) + N" [1 
<g(x)+NS[l 

= q(x) + vl 



l{X t =0} 6' 
1 {Xt=0}] 



J+°°ds X B (<p) 



where we used (|22p for the third equality. This proves ([46 p . 
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Using the Williams' decomposition under N°, we get: 

d r v° r dr N°'W 



VAX) 



J+°°ds X.{ V ) 



Using again the Williams' decomposition under N ,, we have 



J+°° ds X s (tp) 



(50) 



E T 



2/ r d S N°,_ ( c /o + °° dt ^^)-l)l { 



X 3 =0} 



2 f Q ds N° Ys (Jo +oc * x tM -i)l {Xr _ s=0} 



We deduce that, for fixed x, r i— > N 



is non-decreasing and continuous as 



N°[# max = t] = for t > 0. Therefore, we deduce that for fixed x, v t (x) is of class C 1 in t: 



J+°°ds x.( v ) 



d t v t (x) = N°'W 
We have thanks to item (iii) from Proposition 13.71 
(51) No[(e/o +0 °^^(0_i)i {Xs=0} " 

= N° [X s / 0] + N° po + °° d * -1 

= v° s + q{y) - v s {y) 

= \\dx^(v Q s )-d^{y,v s {y)) 



N 



J+ ao dtx t ( v ) 1 



{**#0} 



where the last equality follows from (|15p . (|18p and f)19|) . Thus, with S s (y) = d\ip°(v®) 
d\ip(y,v s (y)), we deduce that: 



N o«) 



J +o ° ^ X.fo,) 



Jo* da s s (y t _ s ) 



This implies (j47j) . Notice that, thanks to (146|) . £ is non- negative and bounded from above by 
2g. □ 

Fix h > 0. We define the probability measures P^ absolutely continuous with respect to P 
and P on T>h with Radon-Nikodym derivative: 

(52) 



Notice this Radon-Nikodym derivative is 1 if the branching mechanism ip is homogeneous. 
We deduce from (@7D and (JMJ) that: 



dP 



(h) 



dF x \V h 



\V± = d h V° h f h dr (8 x $(Y r ,v h - r (Yr))-dx4fl(<<%_ T )) 



dhVh(x) 



and, using f)14[) : 



(53) 



dP 



(h) 

x \T> h 



dPx \v h a ( Y h) d h v h (x) 



9hV h e -J h dr (d^(Y r ,v h _ r (Y r ))-d^°(v°_ r )) 
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In the next Lemma, we give an intrinsic representation of the Radon-Nikodym derivatives 
(|52p and (|53p . which does not involve (3q or v . 

Lemma 4.10. Assume (H1)-(H3). Fix h > 0. The processes = (M t (h) ,t G [0,/t)) and 

MW = (M^ h \t £ [0,/i)), with: 

M (h) _ d h V h - t (Y t ) ft ds dx ^(Y S) v h _ s {Y a )) , jrUh) _ d h v h - t (Y t ) f t ds s ^ ( y sA _ s( y s)) 

d h v h (x) 1 d h v h (x) 

are non-negative bounded T>t -martingales respectively under P x and P x . Furthermore, we have 
for0<t< h: 



(54) 



dP 



(h) 

x \V t 



dP n 



M^ h) and 



dP 



(h) 



'x \V t dP x |D t 

Notice the limit M^ 1 ' of and the limit Ml' 1 ' of M^ 1 ) are respectively given by the 

right-handside of (f53|) and ([52]) . 

Remark 4.11. Comparing ([TO]) and ([34]) . we have that Pjjp = P x with g(t,x) = dhVh-t{x), if 5 
satisfies the assumptions of Remark 13.31 

Proof. First of all, the process is clearly P t -adapted. Using ([4"T|) . we get: 

d h v h - t (y) 



We set: 



We have: 



M 



00 



E, 



1 / h Sfc_,(y.)dr' 



E x [Mf } |A] 



J *E h _ r (y r )dr 



j;s fe _ a (y s )* 



Ey t 



dhVh(x) d h vl_ t 

_ d h v h -t(Y t ) e -tid x j(Y s ,v h _ s (Y B j) ds d h v° h c ti d ^o( v °_.) 
dhVh(x) d hV° h _ t 

In the homogeneous setting, v° simply solves the ordinary differential equation: 

d h v° h = 

This implies that 

d h log(d h v h ) = ^4 = -d^ (v° h ) 

and thus 

(55) ^L e I^^{vl_ B ) d s = 1 

We deduce that 



t/s 



dhvl 



E [Af (h) lP,l - d hVh-t{Y t ) r* drdx $( Yr ,v h _ r (Y r )) 



M t (h) . 
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Therefore, AfC 1 ) is a 2^-martingale under and the second part of (|54|) is a consequence of 
(I52p . Then, use (j!4H to get that is a Pj-martingale under P x and the first part of (|54H . □ 

We now give the Williams' decomposition: the distribution of (Hmax; ^r max ,-Rf' >^T m ) 
under N x or equivalently under N x /a(i;). Recall the distribution P^ defined in (j52[) or ()53|) . 

Theorem 4.12 (Williams' decomposition under N x ). Assume (H1)-(H3). We have: 

(i) T/ie distribution of H max under N x zs characterized by: N x [H max > h] = Vh(x). 

(ii) Conditionally on {H max = ho}, the law ofWx max under N x is distributed as Yt ^ under 
p(M 

(iii) Conditionally on {H max = h^} and WT max , Rj* and Rrp max are under independent 
Poisson point measures on R + x Cl with intensity: 

Mo,h )( s ) ds l{H m «(W")<ho-«} Q! ( w r milB (a)) N w Tmm (s)[dW'}. 
In other words, for any non-negative measurable function F, we have 



F{H max ,W Tmax ,R°& 



max -'max' 



d h v h {x) dh E x h) 



where under E x and conditionally on Yt fA, R w >( h )>9 an d ftW,(h),d are ^ wo i n d e p en d en t Poisson 
point measures with intensity: 

(56) v w >( h \ds,dW) = l m {s)ds l {Hmax (w)<h-sMY s ) N Y3 [dW}. 

Notice that items (ii) and (iii) in the previous Proposition imply the existence of a measurable 
family (N^ , h > 0) of probabilities on (p,,Q) such that is the distribution of W (more 
precisely of (WT max >-Rf> >-^r m )) under N x conditionally on {H max = h}. 

Proof. We keep notations introduced in Proposition 14.71 and Theorem 14.121 We have: 



N, 



F(H m ^W T ^,R 9 T _,R d _] 



N 



e f+°° ds x siv) F(Fmax) W T ^,R^RtJ 
F(H max , W T ^R 9 T _,R d Tm J e^+ R w)(/) 



d h vl dh E x 

oo 

d h vl dh E x 

oo 

d h vl dh E x 

oo 

d h v% dh E x 



F(h, Y [0M , R w >W>s, R w ^ d ) e (R w ^+A w ^' d )(f) 



F(h, Y [0A) ,R W ^, R w >W> d ) e 2/o dS ^ 



=0} 



J- h E fc _.(yi)d« 



F(h,Y l0ih) ,R w ^,R w ^ d ) 



where the first equality comes from (HI) and item (iii) of Proposition 13. 7( we set f(s,W) = 
Jq + °° X r (W)((p) for the second equality; we use Proposition 14.71 for the third equality; we use 
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Lemma 14.31 for the fourth with R W '( h )<9 and R w >(ty> d which under E^ and conditionally on 
Y[o,h) are two independent Poisson point measures with intensity i/W'W; we use (|5ip for the 

fifth, definition (j52j) of Ex for the sixth, and (14"T|) for the seventh. Then use (f33|) and (jT7|) to 
conclude. □ 

The definition of N?' gives in turn sense to the conditional law = N x (.|i/ max = h) of the 
(£, (3, a) superprocess conditioned to die at time h, for all h > 0. The next Corollary is then a 
straightforward consequence of Theorem 14.121 

Corollary 4.13. Assume (H1)-(H3). Let h > 0. Let x € E and Y[o,h) be distributed according 
to P x h ^ . Consider the Poisson point measure M = X^'ejfys-.X.j) on [0, h) X O u>ii/i intensity: 

21 [0jh) (s)ds l{H max {x)<h-s}u{Y s ) n Ya [dX\. 
The process = (X\ h \t > 0), which is denned /or all t > by: 

J, Sj<t 

is distributed according to N^. 

We now give the superprocess counterpart of Theorem 14.121 

Corollary 4.14 (Williams' decomposition under P^). Assume (H1)-(H3). We have the follow- 
ing result. 

(i) Sample a positive number ho according to the law of H max under F u : F u (H max < h) = 
e -^K)_ 

(ii) Conditionally on ho, sample xq 6 E according to the probability measure 

dhVhoix K (dx) 

(iii) Conditionally on ho and xq, sample X^ h °^ according to the probability measure Nio ■ 

(iv) Conditionally on ho, sample X' , independent ofxo and X^ h °\ according to the probability 
measure F u (.\H max <h ). 

Then the measure valued process X' + X^ h °^ has distribution P„. 

In particular the distribution of X'-\-X( h °' conditionally on ho (which is given by (ii)-(iv) from 
Corollary I4.14p is a regular version of the distribution of the (£,/3,a) superprocess conditioned 

to die at a fixed time ho, which we shall write F^ \ 

Proof Let fj, be a finite measure on M + and / a non- negative measurable function defined on FL + x 
E. For a measure-valued process Z = (Z t ,t > 0) on E, we set Z(ffx) = j f(t,x) Z t (dx) fj,(dt) . 
We also write f s (t, x) = f{s + t,x). 

Let X' and X^ be defined as in Corollary 14.141 In order to characterized the distribution 
of the process X' + X^ ho \ we shall compute 



A = E[e- x 'W- xlh0) W]. 
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We shall use notations from Corollary 14,131 We have: 
A = - 



+00 







u(d h v h )e-^dh I ^^l u(dx) 
J E v(d h v h ) 



v{dx 



Y [0,h)] 



+00 



d h v h {x) dh 



E^\E[e-^ xj ^\Y m ] 



■xVti \H <h 



1 {H max <h} 



where we used the definition of X' and TV for the first equality, and the equality ¥ u (H m 
h) = Pj/(-£f m ax < h) = e~ u ( h ^ for the second. Recall notations from Theorem 14.121 We set: 



< 



Vie/ i'ar J 



andg(h)=E v [e~ x (M l {Hn 



c<h}\ 



We have: 



.4 



v(dx) 



+00 



d h v h (x) dhE^ \G(R w ^ 9 ,R w ^' d )g(h) 



is(dx)N x G(R^,R^)g(H max ] 



v(dx) 



\h=H n 



E 



n ' 



{ ^max < ^maxl 



E r e -E ie z^(/M)i 



where we used the definition of G and g for the first and third equalities, Theorem 14.121 for 
the second equality, the master formula for Poisson point measure Yliel ^X* W1 th intensity 
v{dx) ^ x [dX] for the fourth equality (and the obvious notation = inf{i > 0; = 0}) and 

Theorem 12.31 for the last equality. Thus we get: 



This readily implies that the process X' + X^) is distributed as X under V u . 



□ 



5. Some applications 

5.1. The law of the Q-process. Recall F„ defined after Corollary 14.141 is the distribution 
of the (£, /3, a)-superprocess started at v £ A4f(E) conditionally on {H max = h}. We consider 

also Pl- h) = P„( • \H max > h) and N x - h) = N x ( ■ \H max > h) the distributions of the (£,/3,a)- 
superprocess conditionally on {-ff max > h}. 

»(>*) 



The distribution of the Q-process, when it exists, is defined as the weak limit of F u ~ when 
goes to infinity. The next Lemma insures that if 
this limit is also the distribution of the Q-process. 



h goes to infinity. The next Lemma insures that if Fl h) weakly converges to a limit P&> , then 
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Lemma 5.1. Fix t > 0. 7/P^ ) (resp. N ( x h) ) converges weakly to P^, 00 ^ (resp. N x °°^ ) on (Q,J-t), 
then¥^ h ^ (resp. N x ~ h ^ ) converges weakly to p[/°°^ (resp. N x ) on (Q,Ft)- 

Proof. Let Z = 1a with A G such that p[/°°^(cL4) = 0. Using the Williams' decomposition 
under Pj, given by Corollary 14. 14| we have for h > t: 

/•oo 

E^[Z] = e u{vh) / EM[Z] f(h')dh', 

Jh 

where f(h) = — v(dhVh)exp(—i i (vh)). We write down the difference: 

/•oo 

E^ h \z] - E^[Z] = e v{Vh) / (eW[Z] - E^[Z}) f(ti)dti. 

Jh 

Since ^ weakly converges to p[/°°^ on (ft, Ft) and since ¥^°\dA) = 0, we deduce that 

lhn h ,^ +00 El h,) [Z] - E^Z] = 0. We conclude that lim h ^ +00 E [ ^ h) [Z] - E&^Z] = 0, which 
gives the result. The proof is similar for the conditioned excursion measures. □ 

We now address the question of convergence of the family of probability measures 
Recall from flM} that for all < t < h: 

dP {h) 

x l p * = M j h \ 
dPx \v t 

We shall consider the following assumption on the convergence in law of the spine. 

(HA) For all t > 0, P^-a.s. (M^ h \h > t) converges to a limit say M t (oo) , and F, X [M^] = 1. 

Note that Scheffe's lemma implies that the convergence also holds in L 1 (P X ). Furthermore, 

since (M^ h \t G [0, h)) is a non-negative martingale, there exists a version of (M^°°\t > 0) which 
is a non-negative martingale. 

Remark 5.2. We provide in Section [7] sufficient conditions for (Hl)-(HA) to hold in the case of 
the multitype Feller diffusion and the superdiffusion. These conditions are stated in term of the 
generalized eigenvalue Ao defined by 

(57) A = sup {£ G R, 3u G £>(£), u > such that (J3 - C)u = £u}, 

and its associated eigenfunction. 

Remark 5.3. The family (P^ , h > 0) and the family (P x B ' h \ h > 0) defined in Remark 14.61 will 
be shown in Lemma 16.131 to converge to the same limiting probability measure. 

Under (HA), we define the probability measure Pi°°^ on (D,T>) by its Radon Nikodym deriv- 
ative, for all t > 0: 

dp<°°> 



(58) 5p^ = m < 



x \ v t _ ^(oo) 
x \V t 



By construction, the probability measure P x h ^ converges weakly to Pi°°^ on T>t, for all t > 0. 
Let v G A4f(E). We shall consider the following assumption: 
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(H5) u There exists a measurable function p such that the following convergence 
holds in L 1 (z^): 



dhVh 



-> p. 



v(d h v h ) h-^+oo 

In particular, we have v(p) = 1. Let v £ Aif(E). Under (iJ4) and (H5) u , we set: 

p(°°)(dr)= / i/(daOp(x)Pjj.°°)(dy). 



Notice then that j E v(dx) ^^P { x h) (dY) converges weakly to P^°\dY) on V u for all t > 0. 
Remark 5.4. If v a constant times the Dirac mass t^j for some x £ E, then (H5) u holds if (HA) 



holds and in this case we have P 



(oo) 



p1°°^ 



We can now state the result on the convergence of . 

Theorem 5.5. Assume (Hl)-(HA). Let t > 0. The triplet ((W Tm J m , (R 9 Tm J[o,t], i R T m J[o,t]) 
under converges weakly to the distribution of the triplet (Y[ 0jt ], R^'q, R^'i) where Y has 

distribution Pi°°^ and conditionally on Y , R B,d and R B,d are two independent Poisson point 
measures with intensity v B given by 113 9\) . We even have the slightly stronger result. For any 
bounded measurable function F, we have: 



(59) NW 



F ( ( W T ma J [0,t] . ( R T ma J [0,t] » ( R T ma J [0,t] ) 



> 



e(°°) 



Proof. Let h > t. We use notations from Theorems I5.5l and l4.12i Let F be a bounded measurable 
function on W x (M + x f}) 2 . From the Williams' decomposition, Theorem 14. 12} we have: 



F((W Tm J m AK m Jlo,^(^J[o,t]] 



p(V, f> W >9,{h) r>W,d,(h) 



L [o,t] 



where 99™ is defined by: 



/(y [M )=EW 



We also set: 



[o,t] '-"-[o.t] 

F (y[o,t].^»^p;J 



y = y 



^( y[ o, t] )=E(°°) 
We want to control: 



y = y 



Notice that: 



e(°°) 



F ( Y [oAi R [o%' R [o,t}> 



A h = ew [^(y [0 , t] )] - e(°°) [^°°(y[o, t ])] 
(60) = (ew [^(y [0)t] )] - 4 00> + [(/ - O (%*])] • 

We prove the first term of the right hand-side of (|60p converges to 0. We have: 



Ei h) " E i°° ) K(*W])] = E<~> [(M^ - M^) ^(y [( 



r(°°)\ ,Jii 



OA' 
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Then use that (p h is bounded by H-FH^ and the convergence of (Alj: h \h > t) towards Mj 00 "* in 
L 1 ^) to get: 

(61) lim EW [<p h (Y [0A )] - E(°°) K(y [0li ])] = 0. 

We then prove the second term of the right hand-side of (j60|) converges to 0. Conditionally 
on Y, R^'q and R^fj (resp. R^'^ and J^'^) are independent Poisson point measures with 
intensity l[o,t](s) v w '( h \ds,dW) where v w ^ h > is given by (foU]) (resp. u B (ds,dW) where 

v B is given by (139|) ). And we have: 

l[o,t](s) ^^(da,^) = l {ffmax ( W )< ft - s} l[o >t ](s) ^(da,dW0. 
Thanks to (H7J) and ([IBj) . we get that: 

/ 1 {// max (W / )>/i-s}l[o,t](s) ^ B (rfs,dV^) = ds a(y s )N ya [H max > h-s] = ds v h _ s (y s ) < +oo. 

J JO JO 

The proof of the next Lemma is postponed to the end of this Section. 

Lemma 5.6. Let R and R be two Poisson point measures on a Polish space with respective 
intensity v and v. Assume that v(dx) = \A{x)v(dx), where A is measurable and v(A c ) < +oo. 
Then for any bounded measurable function F, we have: 

E[F(R)]-E[F(R)}\<2\\F\\ 00 u(A% 

Using this Lemma with v given by lr 0jt i(s) v B {ds,dW) and A given by {H max (W) < h — s}, 
we deduce that: 



We deduce that: 



(^-O(y[o, t] ) 



E M[(/-^)(F M )] 



< 4IIFI 



ds v h _ s (y s 



< 4 ||i? II e (oo) 

— II Moo x 



ds v h _ s (Y s ) 



Recall that (HI) implies that Vh- S (x) converges to as h goes to infinity. Since v is bounded 
(use (fT71) and (|46l) ). by dominated convergence, we get: 



(62) 



lim E^[( ( p h -^)(Y m )]=0. 



Therefore, we deduce from (1601) that lim 



A/j = 0, which gives ([59 

(oo) 



□ 



We now define a superprocess with spine distribution P 

Definition 5.7. Let v £ Aif(E). Assume P^ 00 ^ is well defined. Let Y be distributed according to 
P^°\ and, conditionally on Y, let M = X^jgJ ^(s 3 -,X>) ^ e a Poisson point measure with intensity: 

21 R+ (s)ds a(Y s )N Ys [dX}. 

Consider the process 

X (oo) = (xi°°\t> 0), which is defined for all t > by: 



X 



(oo) 



jeJ, sj<t 

(i) Let X' independent of X^°°^ and distributed according to P„. Then, we write P^ 00 ^ for 
the distribution of X' + X(°°\ 
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(ii) If v is the Dirac mass at x, we write n1° ' 1 for the distribution of A^°°) 



As a consequence of Theorem 15.51 we get the convergence of P£\ We shall write P x n > when 
v is the Dirac mass at x. 

Corollary 5.8. Under (Hl)-(H4), we have that, for all t > 0: 

(i) The distribution n£ converges weakly to n1°°^ on (Q, J-j). 

(ii) The distribution Px converges weakly to P^°°^ on (O, Ft). 

(iii) Let v G Aif(E). If furthermore {Hh) u holds, then the distribution P^ converges weakly 
to Pl oo) on {Q,T t ). 

Proof. Point (i) is a direct consequence of Theorem 15.51 Definition 15.71 and Proposition 13.121 
Point (ii) is a direct consequence of point (i) , Corollary 14.141 and the weak convergence of 



»0O 



>(<>») 



to Pa; as h goes to infinity. 



According to Corollary 14.141 under p£ . , X is distributed according to X' + X^ where X' 
and X^ are independent, X' is distributed according to P^r h ^ and X^ is distributed according 
to 

idx)^^N^[dX]. 
v{o h v h ) 

Assumption (H5) u implies this distribution converges weakly to: 

v(dx) p(x) N x °°^[dX] 

(because of the convergence of the densities in L 1 ^)) on (0,7^) as h goes to infinity. This and 
the weak convergence of P£ h) to P v as h goes to infinity gives point (iii). 

Proof of Lemma \5.6l Similarly to Lemma 14.31 (formally take / = — ool^c), we have: 



□ 



E [F{R)1 



{R(Ac)=0}\ 



E 



F(R) 



-v{A c ) 



We deduce that: 

E[F(R)] - E[F(R)] | = \e[F(R)] - E[F(R)l {R{Ac)=0} ] e^ 

< \E[F(R)] - E[F(R)1 {R{AC)=0} ]\ + \e[F(R)1 {r{Ac)=0} ](1 - ^ A °> 

< || F \U1 - P(R(A C ) = 0)) + || F IU P(R(A C ) = 0)(e^ AC ) -1) 

= 211^11^(1-6- 
<2\\F\\ 00 v(A c ). 

This gives the result. 



□ 



5.2. Backward from the extinction time. We shall work in this section with the space 
D~ = D(M~,E) equipped with the Skorokhod topology. We also consider the cr-fields T>i = 
o~(Y r ,r 6 /) for I an interval on (— oo,0]. 

Let us denote by 6 the translation operator, which maps any process R to the shifted process 
9 h (R) defined by: 

9 h (R). = R. +h . 
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The process R may be a path, a killed path or a point measure, in which case we set, for 
R = 12jeJ S (s j ,x j ), Oh(R) = Hj&j^ih+Bj^y We also denote P( _h) the push forward probability 
measure of by 6 h , defined on P[_^ )0 ] by: 

(63) P(" ft )(y G •) = ? W (e h (Y) G •) = PW((*W G hMD G •). 

We introduce the following assumptions. 

(HQ) There exists a probability measure on (U _ ,I>(_ DO) o]) denoted p(-°°) such that 
for all x G -E, i > 0, and / bounded and fr_t o] measurable: 

E(- ft ) [/(y ht)0] )] E(-°°)[/(y ht)0] )]. 

(HI) For all i > 0, there exists a non negative function g such that for all x G E, for 
all h > 0: 

/oo 
dr gf(r) < oo. 

Note that the probability measure p( _0 °) in (HQ) does not depend on the starting point x. 

We can now state the result on the convergence of the superprocess backward from the 
extinction time. 

Theorem 5.9. Under {HI) -(HA) and (H6). 

(i) The distribution of the triplet {O h (W Tm J[-t,o], d h(Rg mw )[-t,O],0h(Rd ma *)[-t,o}) under N { x h) 
converges weakly to the distribution of the triplet (Y[_ t0 ]> ^r^o]> rf^to]) w ^ iere Y has 
distribution p(~°°) and conditionally on Y , R w,g and R w,d are two independent Poisson 
point measures with intensity: 

1 {s<o}a(Y s ) ds l {Hmax (W)<-s} N Ys [dW]. 

We even have the slightly stronger result. For any bounded measurable function F , we 
have: 



(64) 



i ? (^(^T mo J [ - t ,0],^(^-)[- t ,0],^(^I'-)h 4 ,0] 



F { Y [-t,0] > R l-t,Q] > ^[4,0] ) 



h— >+oo 

(h) 



(ii) If furthermore (HI) holds, then the process 0/j,PO[-t,o] = (Xh+s,s G [— i,0]) under N x 
weakly converges towards X^^), where for s < 0: 



X (-oo) = ^ .Yi Sj 
is J, Sj<s 

and conditionally on Y with distribution p(~°°) ; XljeJ ^(s-,Xi) * s a Poisson point measure 
with intensity: 

Remark 5.10. We provide in Lemmas 17.31 and 17.61 sufficient conditions for (HQ) and (-HT) to 
hold in the case of the multitype Feller diffusion and the super diffusion. These conditions are 
stated in term of the generalized eigenvalue Aq defined in (|5"T|) and its associated eigenfunction. 
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Proof. Let < t < h. We use notations from Theorems 14.121 15-51 and 15,91 Let F be a bounded 
measurable function on W~ x (R _ x fl) 2 with W~ the set of killed paths indexed by negative 
times. We want to control 6h defined by: 



6 h = N<P 



^(^(WT maj£ )[- t) o],^(i?La,)H,o],^(^ maj£ )H,o]) 



Wc set: 



T(y Hi0] ) = E(—) 



* \.y[-t,o]> H[- t ,o]> K [-t,o] 



e(~°°) 

) 



F(Y 



l-t,0], -K[- t! o]' U [-t,0]) 



Y = y 



We deduce from Williams' decomposition, Theorem I4.12| and the definition of R w,a and R w,d , 
that: 



F{e h (w Tm j hm , o h (K m Ji-t,o], o h (i% m J l - m ) = E(- ft ) [T(Y hm )] . 



We thus can rewrite 5h as: 

S h = E^[T(Y hm )] -E(-°°)[T(y H! o])]- 

The function T being bounded by ||-F||oo an d measurable, we may conclude under assumption 
(H6) that lim/j_5. +OC) 5h = 0. This proves point (i). 

We now prove point (ii). Let t > and e > be fixed. Let F be a bounded measurable 
function on the space of continuous measure- valued applications indexed by negative times. For 
a point measure on M~ x 0, M = Yliex ^(st.Wj); we se ^ : 

^ = ^((£M*W)))hJ- 



For h > t, we want a control of 5^ defined by: 



4 = Ni h ) 
By Corollary I4.13( we have: 



F{e h (x) { _ m ) 



E(-o°) 



Thus, we get: 

(65) ^ = 



-t,0], 



e<-°°) 



For a > s fixed, we introduce 5%, for h > a, defined by: 



(66) 



St = N<*> 



F(^(^ max + < na j ha , 0] ; 



e(-°°) 



F(R W '<> + R w > d ) 



F{(R W * + R w ' d ) { _ at0] ) 



Notice the restriction of the point measures to [—a, 0]. Point (i) directly yields that lim/ l _ > _|_ 00 5% 
0. Thus, there exists h a > such that for all h > h a , 

5 a h < e/2. 
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We now consider the difference 5h — Si- We associate to the point measures M introduced above 
the most recent common ancestor of the population alive at time —t: 

A{M) = sup{s > 0; 1 {s,<~s}l{H n ^(w i )>-t~s l } ¥= 0}. 
iex 

Let us observe that: 

(67) NW a.s., F(9 h (I^ w + R d Tm J)l {A < a} = F{9 h {R 9 T _ + < nax )[-a,o])l { A<a } , 

with A = A(6h(Rj, + Rt^ )[—h,0]) i n the an d i n the right hand side. Similarly, we have: 

(68) p(-°°) a.s., F[R W ' 9 + R W ' d )l {A <a} = F{(R W ' 9 + R W4 )[-a,o])l { A<a}, 

with A = A{R w,g + R W ' ) in the left and in the right hand side. We thus deduce the following 
bound on 6h — S%: 



\S h 



%\<2\\F\ 



= 2 \\ F \ 
< 8\\F\ 



NjW [A > a] + p(-°°) [A > a 

E^T^ [l — e~ fa dr 2Q K-t--"r)(y_ r )j _|_ g(-oo) j-j _ e -f™dr 2a(v r _ t —u r )(Y'_ r )j 
dr g(r), 



a—t 



where we used (|65j) . (|66j) . (j67j) and (j68j) for the first inequality, the definition of A for the first 
equality, as well as {H7) and the fact that 1 — e~ x < x if x > for the last inequality. From (H7), 
we can choose a large enough such that: \5h — 8%\ < e/2. We deduce that for all h > max(a, h a ): 
\8~h\ ^ \$h — 5~h \ + \^h\ — £ - This proves point (ii). □ 

6. The assumptions {HA), {Hh) v and {HQ) 

We assume in all this section that P is the distribution of a diffusion in ~R K for K 
integer or the law of a finite state space Markov Chain, see Section [7] and the references 
therein. In particular, the generalized eigenvalue Ao of {(3 — C) (see (I86p or (I88p ) is known to 
exist. We will denote by 4>q the associated right eigenvector. We shall consider the assumption: 

{H8) There exist two positive constants C\ and C2 such that Vx EE, C\ < (po{x) < C2; 
and <p Q G V{£). 

Under {H8), let Pp be the probability measure on {D,T>) defined by ([9]) with g replaced by 

4>o- 



(69) 



V* > 0, 



x \V t 



dP 



x \V t 



MYo) 



We shall also consider the assumption: 
{H9) The probability measure P^ admits a stationary measure it, and we have: 



(70) 



sup 

/£6£,ll/||oo<l 



°[f{Yt)]-n{f)\ 



t—t+oo 



-> 0. 



Notice the two hypotheses {H8) and {H9) hold for the examples of Section [71 see Lemmas 
1731 and 1731 



Let us mention at this point that we will check that P 



00 



Pl oo) with Pi oo) defined by (I55]>. 



see Proposition [67 
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6.1. Proof of (H4)-(H6). Notice (H9) implies that the probability measure Pt° admits a 
stationary version on D(M., E), which we still denote by P^°. 

We introduce a specific /i-transform of the superprocess. From Proposition 13.51 and the defini- 
tion of the generalized eigenvalue (|86p and (j88[) , we have that the /i- transform given by Definition 
13.41 with g = (f>o of the (£,(3, a) superprocess is the (C^°, Ao, a4>o) superprocess. We define 
for all t > and x G E by: 

vf°(x) = Nl c *°> Xo ' aM \H„ 



(71) 



>t]. 



Observe that, as in (|17p . the following normalization holds between and 



v: 



(72) 



«*»(*) 



vt(x) 
4> (x) 



Our first task is to give precise bounds on the decay of vf° as t goes to oo. 

We first offer bounds for the case Ao = in Lemma [6. H relying on a coupling argument. This 
in turn gives sufficient condition under which (HI) holds in Lemma [6. 21 We then give Feynman- 
Kac representation formulae, Lemma 16.31 which yield exponential bounds in the case Ao > 0, 
see Lemma 16.41 We finally strengthen in Lemma 16.61 the bound of Lemma 16.41 by proving the 
exponential behavior of vf° in the case Ao > 0. The proofs of Lemmas 16.11 16.21 16. 3| 16.41 and 16.61 
are given in Section 16.21 

We first give a bound in the case Ao = 0. 

Lemma 6.1. Assume Xq = 0, (H2) and (H8). Then for all t > 0: 



a<po(x) 



lac 



2 

oo 



< t vf°(x) < a4> (x) 



O 



A coupling argument then implies that (HI) holds: 
Lemma 6.2. Assume A > 0, (H2) and (H8). Then (HI) holds. 

We give a Feynman-Kac's formula for Vq and dv^. 
Lemma 6.3. Assume A > 0, (H2)-(H3) and (H8). Let e > 0. We have: 



(73) 
(74) 



J h+e 



,<k) 



(x) 



e -\ h E 0o 
e -X h E O 



- Jo" ds a(Y s ) MY*) ^+ £ - s ( y s) v fo (y h ) 
-2 J Q h ds a(Y s ) MYs) vtl e _ s (Y s ) g^fo ^ 



We give exponential bounds for Vq and dtv^ in the subcritical case. 

Lemma 6.4. Assume Ao > 0, (H2)-(H3) and (H8). Fix to > 0. There exists C3 and C4 two 
positive constants such that, for all x G E , t > to: 

(75) C 3 < vf°(x)e Xot < C 4 . 

There exists C5 and Cq two positive constants such that, for all x G E, t > to-' 

(76) C 5 <|^f°(x)|e Aoi <C 6 . 

As a direct consequence of ()75|) . we get the following Lemma. 



Lemma 6.5. Assume A > ; (H2)-(H3) and (H8). Then (H7) holds. 
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In what follows, the notation 0^(1) refers to any function such that lim/j^+oo || || = 0. 
We now improve on Lemma 16, 4] by using the ergodic formula (|70p . 



Lemma 6.6. Assume Xq > 0, (H2)-(H3) and (H8)-(H9) hold. Then for all e > 0, we have: 

(77) d t vfl £ (x) e Xoh = Et°[e~ 2 ^ ds a *° v t+e<y->) d t vp(Y )] (1 + o h (l)). 
In addition, for Xq > 0, we have that: 

Et° [ e" 2 Jo°° ds a ^ vtU Y -°) d t vp (F )] 
is finite (notice the integration is up to +oo) and: 

(78) d t vt° +e (x) e^ h = Et[ e - 2 ^ dsa ^ v "U Y -^d t vp{Y Q )]+o h {\). 



Our next goal is to prove {HA) from (H8)-(H9), see Proposition 16.81 

Fix x G E. We observe from (|54p and (j69[) that P^ is absolutely continuous with respect to 
Pt° on ^[o,t] for < t < h. We define the corresponding Radon-Nikodym derivative: 



ciP (/l) 

M (h),<h = x l p [Q,t] 

1 P[o,t] 

Using (|54p . (|69p and the normalization v(x) = v^°(x) 4>o(x), we get: 

(7q) M (hUo = d t V h . t (Y t )e- X ^ MYo) -2f*ds a(Y s ) v h _ a (Ys) 

{ ' 1 " dtv h (Y ) MYt) 

= dtv^jYt)^ 1 ^_ 2 jt ds a{Ys) MYs) v , _ AYs) 
dtvt (Y ) 

We have the following result on the convergence 
Lemma 6.7. Assume (H2)-(H3) and (H8)-(H9). For A > 0, we have: 

M (h),4>o y 1 p^o. as and inL 1 ^ ), 

h—t+oc 

and for Xq > 0, we have: 



M (hUo y i p<f>°-a.s. andinV-CP*? 



Proof. We compute: 



* d^ (Y )e^ 



T ( 1 + ( 1 ) ) 



l + o h (l), 



where we used ([79]) for the first equality, ([77]) twice and the boundedness of a and 0o as well as 
the convergence of Vh to for the second, and Lemma 16.41 (if Xq > 0) or Lemma 16.11 (if Aq = 0) 
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for the fourth. Since Oh(l) is bounded and converges uniformly to 0, we get that the convergence 
of M\ h towards 1 holds P* -a.s. and in L l (?t )- 

Similar arguments relying on ([751) instead of ([77]) imply that M^^° = 1 + 0^(1) for Ao > 0. 

Since o^(l) is bounded and converges uniformly to 0, we get that the convergence of M^fy^ 
towards 1 holds p£°-a.s. and in L^P* ). □ 
The previous Lemma enables us to conclude about (Hi). 



(oo) 

X 



Proposition 6.8. Assume A > 0, (H2)-(H3) and (H8)-(H9). Then (HA) holds with P 
Pp. 

Proof. Notice that: 

M (h) = * \V [0A = M (h),j>o * |P[o,ti 
dP x P[o,t] dF * [D[o,H 

The convergence limh_;. +00 M^'^° = 1 P^°-a.s. and in L 1 (Pf°) readily implies (HA). Then, use 
USED to get p(°°) = P^°. □ 

Notice that (H5) u is a direct consequence of Lemma 16.61 

Corollary 6.9. Assume A > 0, (H2)-(H3) and (H8)-(H9). Then (Hh) v holds with p = 
4>o/v((t>o)- 

Proof. We deduce from J7J) and (JTZJ) that: 

^(x) = /(fc)0o(a?) (1 + o h (l)) e~ Xoh , 
for some positive function / of h. Then we get: 

dhVh(x) M x ) „ , „ /-m 
v[d h v h ) v(<po) 

This gives (H5) u , as 0^(1) is bounded, with p = (j>o/v(4>o)- D 
Our next goal is to prove (HQ) from (H8)-(H9), see Proposition 16. 121 

Observe from (|53p . ()63p and (|69p that P^- ^ is absolutely continuous with respect to P^° on 
Pr-^oj- We define L^~ h ^ the corresponding Radon- Nikodym derivative: 

(80) L(~V = |P[ - fc ' 01 = 1 dhV h c ~2 r° h (a(y J )t>- J (Y«)-«° J ,) da 

dP*° m a(y )<fo(*b) ^"(^Oe^ 

The next Lemma insures the convergence of L^~ h ^ to a limit, say L(-°°). 
Lemma 6.10. Assume A > 0, (H2)-(H3) and (H8)-(H9). We have: 

L i-h) y L (-oo) p0o. fl S a^fflL^pJo). 

ft— >+oo 

Proof. Notice that lim^ +t>0 d h v° h = -ffi. We also deduce from (gSJ), (09]) and {75]) that 
f® h (a(Y s )v- s (Y s ) — v^_ s ) ds increases, as h goes to infinity; to (a(Y s )v- s (Y s ) — v°_ s ) ds which 
is finite. For fixed t > 0, we also deduce from (|78p (with h replaced by h — t and e by t) that 
P^r a.s.: 

lim a^°(y_ /l )e Aoh = e Xot Et°\e- 2 ^ ds a{Ya) * o( - Y ^ v -°^ d t vf° (Y- t )] . 

ft.— >+oo 
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We deduce from ([80]) the a.s. convergence of (L(~ h \h > 0) to L0-°°). Notice from ||7Bjl 
that, for fixed t, the sequence {L^~ h \h > i) is bounded. Hence the previous convergence holds 
also in L^P* ). □ 

As Et° [lHO] = 1, we deduce that e£° [L^ 00 )] = 1. We define the probability measure 
P^r co ^'^° on (D~ ,T>r OQ rt) by its Radon Nikodym derivative: 

^p(-°°).0o 

(81) 1 P '-' 01 =L(-°°). 

dP^% 

T l^(-oo,0] 

Remark 6.11. Assume A > 0, {H2)-{m) and (#8)-(#9). Define for ft > t > 0: 

dP ( ~ h) 

L^) = E^[L(^)|P ( _,_ t] ] = -^^ 

«P IT) 

,p(-Oo),0 O 

L (_-oo) = Ejfo^-oo)^^] = ^ |P <--*' ■ 

T P(-oo,-t] 



Using ([55]) and Lemma [673l we get: 

= a^t(y-t) ^o(^-fe) -Ao(fe-t) -2 Cfeds a(y s )«_ s (y s ) 
-* $i; h (Y_ h ) fo(Y_ t ) 

" [ e -2 ^ ds a(Ys) ^ o(ys) ^ (ys) 9ivf° (YL t )] 

Using Lemma [6.61 and convergence of (L { ~ t h \ h > t) to L^* , which is a consequence of Lemma 
16.101 we also get that for t > 0: 

^ ( _ oo) _ e -2,Q g 0o(^)> s (^)^fo ( y_ t) 

"* ~ Et° [ e" 2 ds a{Ys) MYs) v -^ Ys) d t vf° (F_ t )] 

Those formulas are more self-contained than ([80[) and the definition of L^~°°^ as a limit, but 
they only hold for t > 0. 

The following Proposition gives that {HQ) holds. 

Proposition 6.12. Assume A > 0, (H2)-(H3) and (H8)-(H9). Then (H6) holds withP^ 00 ^ = 

p(-oo),0 o _ 

Proof. Let < t and F be a bounded and V \-t y o] measurable function. For h large enough, we 
have: 

4" fc) [F(Y { _ m )} = E W [E^f [F(^ /2 (y) H ,_ s] )]] 

V 2 



E*[ilf$* E^[ J P(^ /2 (Y) [ _ ti0] )]] 

E^[E^[F(0 V2 (y) HiO] )]]+ o ,(i) 
4 fc/2) [^ fc /2CO[-t,ti|)]+ofc(i) 

E (-VV[F(Y hm )]+o h (l), 
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where we used the definition of p( - ' 1 ) and the Markov property for the first equality, Lemma 
6.71 together with F bounded by ||-F||oo for the third, and assumption (H9) for the fourth. We 
continue the computations as follows: 

E(-V[F(Y { _ m )] =Et°[L(- h MF(Y { _ m )] + o h (l) 
= E £o[ L (-oo ) F{ Y { _ m )]+o h {l) 

= Ei—)^[F(Y [ _ ti o ] )]+ 0/l (l), 

where we used Lemma [6. 101 for the second equality. This gives (HQ) with p(-°°) = p(-°°).*>. □ 

6.2. Proof of Lemmas mi 131 [63], [63] and S 

Proof of Lemma \6.1[ From (#2) and (H8), there exist m, M S M such that 

\/x e E, < m < a<j) (x) < M < oo. 

M_ 

M 



Let W be a (^^£,0,M) Brownian snake and define the time change <E> for every KiGWby 



&t(w) = J ds -M^(w(s)). As dt&t(w) > 1, we have that t — > &t(w) is strictly increasing. Let 

t — > &\ ( w ) denote its inverse. Then, using Proposition 12 of [ID], first step of the proof, we 
have that the time changed snake W o $ , with value 

(W o $-!) s = (W S (^\W S )), t E [0, *-\W a , H,)]) 

at time s, is a (£,0,a</>o) Brownian snake. Noting the obvious bound on the time change 
&^ 1 (w) < t, we have, according to Theorem 14 of jlOj : 

^ ^pjx) , l-Hmax S t) > (,-Hmax S t) 

u Ox 



which implies: 

a 



M 

from the exponential formula for Poisson point measures. Now, the left hand side of this in- 
equality can be computed explicitly: 

(J"_ £<#>() ,o,m) (4r-£ 0o ,o,M) i 

' (#max > t) = ^ >(H max >t) = — 

and the right hand side of this inequality is vf°(x) from ([71]) . We thus have proved that: 

and this yields the first part of the inequality of Lemma [6. 11 The second part is obtained in the 
same way using the coupling with the (-^j^£^°, 0, m) Brownian snake. □ 

Proof of lemma \6.2[ Assumption (H2) and (H8) allow us to apply Lemma 16.11 for the case 
Ao = 0, which yields that vto = 0, and then Voo = thanks to (|72p . This in turn implies that 
(HI) holds in the case Ao = according to Lemma 12.51 For Ao > 0, we may use item 5 of 
Proposition 13 of [ID] (which itself relies on a Girsanov theorem) with p(£'°> a 0°) in the role of P c 
and f>(^ °>*o> a 0o) in the role of ¥ b,c to conclude that the extinction property (HI) holds under 

p(£* ,Ao,a^o). □ 
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Proof of Lemma \6.3[ Let e > 0. The function is known to solve the following mild form of 
the Laplace equation, see equation ([3]): 

f-t 



J t+S 



(x)+Ef° 



dr (Ac vf° s _ r (Y r ) + a{Y r )cf> {Y r ){vil s _ r {Y r )f 



Et°[vf°(Y t )]. 



By differentiating with respect to s and taking t = t — s, we deduce from dominated convergence 
and the bounds ([46]) . ([47]) and ([49]) on = v/<po and its time derivative (valid under the 
assumptions (H1)-(H3)) the following mild form on the time derivative dtv^°: 

ft- 8 



d t vt°(x) + Ei° 



dr (A +2a(Y r )MYr)vt r (Y r ))d t vt r (Y r ) 



Ep[d t vf°(Y t _ s )). 



From the Markov property, for fixed t > 0, the two following processes: 



vf° s (Y s ) - I dr (A + a(Y r )MYr)vtMYr)) vt r (Y r ), 0<s<t 



and 



9tvt s (Y s ) - j* dr (A +2a(Y r )MYr)vt\(Y r ))d t vl° r (Y r ),0<s<t 



are 2? s -martingale under Pf° ■ A Feynman-Kac manipulation, as done in the proof of Lemma 
13.1] enables us to conclude that for fixed t > 0: 



0<s<t 



and 



d t V?° s (Y s ) e" & dT ( A ° +MYr)MYr)vf° r (Yr)) 



0< s <t 



are D s -martingale under P%° . Taking expectations at time s = and s = h with t = h + e, we 
get the representations formulae stated in the Lemma: 



-Ao/i p^o 



□ 



Proof of Lemma \6.4\ Since vf° = v e /4>o = v e /(a(j>o), we can conclude from (|46p . (#2) and (i?8) 
that vf° is bounded from above and from below by positive constants. Similarly, we also get 
from (jlTj) . ([18]) and ([M]) that |^t; E | is bounded from above and from below by two positive 
constants. Thus, we have the existence of four positive constants, D±, D 2 , D3 and D4, such 
that, for all x G E: 

(82) D 1 <v*°(x)<D 2 , 

(83) D 3 <\d t vp(x)\<D 4 . 

From equations ([73]) . ([82]) and the positivity of -y^ , we deduce that: 

(84) vf^ix) <D 2 e~ Xoh . 

Putting back ([841) into ([73]) . we have the converse inequality D§e~ x ° h < vf^ £ (x) with = 
L>iexp{-D 2 Halloo H^ML/Ao} > 0. This gives ([75]). 

Similar arguments using ([71]) and ([83]) instead of ([73]) and ([82]) . gives ([76]) . □ 
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Proof of Lemma \6.6l Using the Feynman-Kac representation of dhvf^ +£ from (|73p as well as the 
Markov property, we have: 



d h vil £ {x)^ h = ^t° 

Notice that 
(85) 



e -2 j* ds a(Y s ) MYs) 4° +£ _ 3 (Y a ) ^ 
e -2 ffds a 0o vi° +s _ s (Y s ) E p [ e -2 j^cfe a 0o <° ^ +£ _ s (^) ^ (y^)] 







< Hal 



OO II ' Hoo 



according to Lemma 16.41 if Ao > and Lemma 16. II if Ao = 0. We get: 



r?d>o 

J-I 7r 
■p</>o 

7T 
E 0O 



'/<() 



+ ^ Js; ^°(y h _^) 



(l + o h (l)) 



where we used ()85[) for the first equality, (#9) for the second, stationarity of Y under Pi for 
the third and ([85]) again for the last. This gives ([771) . 
Moreover, if Aq > 0, we get that: 



- 2 /_« ds a ^° s ( y ») (y ) 



is finite and that: 
lim 

h'—^+oo 



Therefore, we deduce ([75]) from ((7 



□ 

6.3. About the Bismut spine. Choosing uniformly an individual at random at height t under 
N-r and letting t — > oo, we will see that the law of the ancestral lineage should converge in some 
sense to the law of the oldest ancestral lineage which itself converges to Pl oo) defined in (p]h 
according to Lemma 16.81 

We have defined in (|45p the following family of probability measure indexed by t > 0: 



dP 



(B,t) 
x \V t 



-P dsP{Y a ) 



dP 



x \v t E x 



Sods 0(Y S ) 



Lemma 6.13. Assume (H8)-(H9). We have, for every < to < t: 



dP 



(B,t) 
x \V tQ 



dP 



(oo) 

P x -a.s. and in L 1 (P X ) 



dP x \v tQ *^+°° dP x \ VtQ 
Note that there is no restriction on the sign of Aq for this Lemma to hold. 
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Remark 6.14. This result correspond to the so called globular phase in the random polymers 
literature (see [8], Theorem 8.3). 

Proof. We have: 



dP 



B,t) 
x \V tQ 



dP 



x \V tQ 





e - /*-*» ds P(Y S ) 


E x 









e -/ * to ds(/3(y a )-Ao)" 


E x 







- / 4 ° (/3(y s )-A ) M Y to) EYt ° 



e -/ *-*° ds (/3(y s )-A ) Mg^fr) 1 



dP 



E 



00 



[1/MYt-to)] 



E , 



„-/o ds (/3(y s )-A ) 0o (yp 1 
e 0o (x) MYt) 



dP„ 



x \v tQ Et° [i/MM 

t-foo dP x PtQ 7r(i) ~ dP x {VtQ ' 

where we use the Markov property at the first equality, we force the apparition of Ao at the 
second equality and we force the apparition of 4>q at the third equality in order to obtain the 
Radon Nikodym derivative of P^ with respect to P x : this observation gives the fourth equality. 
The ergodic assumption (H9) ensures the P x -a.s. convergence to 1 of the fraction in the fourth 
equality as t goes to oo. Since 

\t,y) [1/MYt-to)) /Ef° [I/MY)] 

is bounded according to (H8), we conclude that the convergence also holds in L 1 (P X ). Then use 
Lemma ESI to get that Pt° = P^ . □ 



7. TWO EXAMPLES 

In this section, we specialize the results of the previous sections to the case of the multitype 
Feller process and of the super diffusion. 

7.1. The multitype Feller diffusion. The multitype Feller diffusion is the superprocess with 
finite state space: E = {1, . . . , K} for K integer. In this case, the spatial motion is a pure 
jump Markov process, which will be assumed irreducible. Its generator £ is a square matrix 
{lij)i<i,j<K of size K with lines summing up to 0, where qij gives the transition rate from i to 
j for i 7^ j. The functions (3 and a defining the branching mechanism ([2]) are vectors of size K: 
this implies that (H2) and (H3) automatically hold. For more details about the construction of 
finite state space superprocess, we refer to [H], example 2, p. 10, and to [6] for investigation of 
the Q-process. 

The generalized eigenvalue Ao is defined by: 

(86) A = sup {£ e M, 3u > such that (Diag(/3) - C)u = £u}, 

where Diag(/3) is the diagonal K x K matrix with diagonal coefficients derived from the vector 
(3. We stress that the generalized eigenvalue is also the Perron Frobenius eigenvalue, i.e. the 
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eigenvalue with the maximum real part, which is real by Perron Frobenius theorem, see [30| . 
Exercise 2.11. Moreover, the associated eigenspace is one-dimensional. We will denote by (fro 
and <f>o its generating left, resp. right, eigenvectors, normalized so that YliLi 4 > o(i)4'o{i) = 1> an d 
the coordinates of 4>o and cfto are positive. 

We first check that the two assumptions we made in Section 6 are satisfied. 

Lemma 7.1. Assumptions (H8) and (H9) hold with tt = <pQ (j>Q. 

Proof. Assumption (H8) is obvious in the finite state space setting. Assumption (H9) is a 
classical statement about irreducible finite state space Markov Chains. □ 

Lemma 7.2. Assume A > 0. Then {HI), (HA) and (Hh) v hold. 

Proof. Assumption (H2) and (H8) hold according to Lemma 17.11 Together with Ao > 0, this 
allows us to apply Lemma 16.21 to obtain (HI). Then use Proposition 16.81 to get (H4) and 
Corollary E3S to get (Hh) u . □ 

Lemma 7.3. Assume Ao > 0. Then (H6) and (H7) holds. 

Proof. We apply Proposition 16.121 to prove (HQ) and Lemma 16.51 to prove (H7). □ 
Recall that and p1°°^ were defined in ([34"]) and ([35]) respectively. 



Lemma 7.4. We have: 

(i) P^ is a continuous time inhomogeneous Markov chain on [0, h) issued from x with 
transition rates from i to j, i ^ j, equal to |jj Qij a ^ time t, < t < h. 

(ii) Pi°°^ is a continuous time homogeneous Markov chain on [0, oo) issued from x with 
transition rates from i to j, i ^ j, equal to ^jjjQij- 

Proof. The first item is a consequence of a small adaptation of Lemma 13.21 for time dependent 
function. Namely, let gt{x) be a time dependent function. Consider the law of process (t,Yt) 
and consider the probability measure P g defined by Q with g(t,Y t ) = gt(Y t ). Denoting by C\ 
the generator of (the inhomogeneous Markov process) Yt under P 9 , we have that: 

(87) VuGP 9 (£), Cl{u) = ^9tu)-C{g t )u 

9t 

Recall that for all vector u, C(u)(i) = Ylj^iQij{ u (j) ~ U W)- Then apply ([87]) to the time 
dependent function gt(x) = dtVh-t(x), and note that P 9 = P^^ thanks to (|54p . For the second 
item, observe that Proposition l6.8l identifies p1°°^ with P^°. Use then Lemma f3.2l to conclude. □ 

William's decomposition under (Propositions I4.14p together with the convergence of 
this decomposition (Theorem 15. 5p then hold under the assumption Ao > 0. Convergence of 
the distribution of the superprocess near its extinction time under (Proposition 15. 9p holds 
under the stronger assumption Ao > 0. We were unable to derive an easier formula for p(~°°) 
in this context. 

Remark that Lemma 15.11 Definition 15.71 and Corollary 15.81 give a precise meaning to the 
"interactive immigration" suggested by the authors in Remark 2.8. of [6]. 
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7.2. The superdiffusion. The superprocess associated to a diffusion is called superdiffusion. 
We first define the diffusion and the relevant quantities associated to it, and take for that the 
general setup from [27]. Here E is an arbitrary domain of M. K for K integer. Let and bi 
be in C l,IJ, (E), the usual Holder space of order \x G [0,1), which consists of functions whose 
first order derivatives are locally Holder continuous with exponent fx, for each i,j in {1, . . . , K}. 
Moreover, assume that the functions a^j are such that the matrix (oi i3 - )(i,j)e{i...K'} 2 i s positive 
definite. Define now the generator C of the diffusion to be the elliptic operator: 

K K 
i=l i,j'=l 

The generalized eigenvalue Ao of the operator /3 — C is defined by: 

(88) A = sup {£ eR,3ue £>(£), u > such that (/3 - C)u = I u}- 

Denoting E the expectation operator associated to the process with generator C, we recall an 
equivalent probabilistic definition of the generalized eigenvalue Aq: 




for any x G R , where ta c = mf {t > : Y (t) ^ A} and the supremum runs over the compactly 
embedded subsets A of R . We assume that the operator (/3 — Ao) — C is critical in the sense 
that the space of positive harmonic functions for (/3 — Ao) — C is one dimensional, generated by 
(fro. In that case, the space of positive harmonic functions of the adjoint of (/3 — Ao) — C is also 
one dimensional, and we denote by (fro a generator of this space. We further assume that the 
operator (J3 — Ao) — C is product-critical, i.e. j E dx <fro(x) (fro(x) < oo, in which case we can 
normalize the eigenvectors in such a way that f E dx (fro(x) (fro(x) = 1. This assumption (already 
appearing in [15]) is a rather strong one and implies in particular that P^ is the law recurrent 
Markov process, see Lemma 17.51 below. 

Concerning the branching mechanism, we will assume, in addition to the conditions stated in 
section El that a £ C 4 (E). 

Lemma 7.5. Assume (H8). Assumption (H9) holds with ir(dx) = (fro(x) (fro(x) dx. 

Proof. We assume that (/? — Ao) — C is a critical operator which is product critical. Note that 
— is the (usual) /i-transform of the operator (/3 — Ao) — C with h = (fro, where the /i-transform 
of £(•) is C{h-)/h. Then Remark 5 of [15] implies that -C*° is again a critical operator which 
is also product critical with corresponding (fro and (fro given by 1 and (fro (fro- Then Theorem 9.9 
p. 192 of [27], see (9.14), states that (H9) holds. □ 

Note that the non negativity of the generalized eigenvalue of the operator (/3 — C) now charac- 
terizes in general the local extinction property (the superprocess X suffers local extinction if its 
restrictions to compact domains of E suffers global extinction); see [16] for more details on this 
topic. However, under the boundedness assumption we just made on a and (fro, the extinction 
property (HI) holds, as will be proved (among other things) in the following Lemma. 

Lemma 7.6. Assume Ao > and (H8). Then (H1)-(H4) and (H5) u hold. If moreover Ao > 0, 
then (H6) and (H7) holds. 

Proof. The assumption a £ C 4 (E) ensures that (H2) and (H3) hold. Then the end of the proof 
is similar to the end of the proof of Lemma 17.21 and the proof of Lemma 17.31 □ 

Recall that P { x h) and p1°° ) were defined in flM! and ([58D- 
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Lemma 7.7. We have: 

• P x ^ is an inhomogeneous diffusion on [0, h) issued from x with generator at time t £ 



• Pi° is an homogeneous diffusion on [0, oo) issued from x with generator iC + a-^V. 



William's decomposition under N x (Propositions Bmj) together with the convergence of this 
decomposition (Theorem [53]) then hold under the assumption Ao > and (H8). Convergence of 
the distribution of the superprocess near its extinction time under (Proposition 15 . 9[) holds 
under the stronger assumption Ao > 0. 

Remark 7.8. Englander and Pinsky offer in [T7] a decomposition of supercritical non-homoge- 
neous super diffusion using immigration on the backbone formed by the prolific individuals (as 
denominated further in Bertoin, Fontbona and Martinez [1]). It is interesting to note that the 
generator of the backbone is C w where w formally satisfies the evolution equation Cw = ip(w), 
whereas the generator of the spine of the Q process investigated in Theorem 15.51 is where 
(f)Q formally satisfies £(fro = (3 4>q. In particular, we notice that the generator of the backbone C w 
depends on both /3 and a and that the generator of our spine does not depend on a. 
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